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THERMAL BUCKLING ANALYSIS FOR STIFFENED 


ORTHOTROPIC CYLINDRICAL SHELLS* 

By L. K. Chang** and Michael F. Card 
Langley Research Center 

SUMMARY 

A theory for thermal buckling of an orthotropic, multilayered, stiffened cylindrical 
shell is presented. The theory includes the effects of eccentricity of layers and stiff- 
ening, and deformations prior to buckling. It is sufficiently general to account for dis- 
crete rings and averaged properties of longitudinal stiffening, as well as arbitrary tem- 
perature distributions through the thickness of the shell and depth of the stiffeners. Two 
computer programs are described corresponding to solutions for buckling obtained by 
using finite differences and determinant plotting or modal iteration. 

Computed results for thermal buckling of unstiffened and ring-stiffened shells are 
presented and are in reasonable agreement with published results. The interaction of 
thermal loading and axial compression in two large-diameter stiffened shells representa- 
tive of a launch vehicle interstage and a preliminary supersonic transport fuselage design 
is investigated. Results indicate that buckling can occur in both structures at a realistic 
temperature under thermal loading alone. 

INTRODUCTION 

Advanced aerospacecraft of the seventies, such as space shuttles and hypersonic 
aircraft, will probably experience severe thermal environments during flight. A renewal 
of interest in thermal stress and thermal buckling problems is anticipated, since it is 
unlikely that thermal protection can completely eliminate stresses and deformations 
arising from thermal gradients in stiffening frameworks, bulkheads, and covers. Ther- 
mal buckling problems in unstiffened shells have been investigated in references 1 to 8. 
Results for stiffened and sandwich shells are given in reference 3 and references 9 to 12. 
An excellent summary of existing thermal buckling analyses and experiments is given in 
reference 13. A review of these references, however, indicates that further studies of 
thermal buckling in stiffened, multilayered shells are needed. 

*Part of the information herein was presented at the AIAA/ASME 11th Structures, 
Structural Dynamics, and Materials Conference. 

**NRC-NASA Resident Research Associate, now at Argonne National Laboratory, 
Argonne, 111. 



Recent advances in computer technology together with derivations of more consis- 
tent buckling theories (refs. 14 to 17) permit the development of buckling analyses which 
accurately account for thermal stresses and deformations in stiffened shells. The pur- 
pose of the present paper is to present an analysis for thermal buckling of stiffened 
cylindrical shells, to describe two computer programs based on solutions for buckling, 
and to report the results of thermal buckling calculations for unstiffened and stiffened 
cylindrical shells. 

The analysis was developed by extending the theory of reference 17. The theory is 
sufficiently general to account for the averaged properties of longitudinal stiffening and 
multiple orthotropic layers with differing coefficients of expansion as well as arbitrary 
temperature gradients through the thickness of the shell and depth of the stiffeners. A 
general outline of the theory is given and buckling solutions suitable for rapid computer 
calculations are presented. Two computer programs based on determinant plotting and 
modal iteration are described. Results of computations for unstiffened and ring-stiffened 
cylindrical shells under thermal loads are compared with results from existing buckling 
theory. Results are presented to demonstrate the effects of interactions of axial com- 
pressive loads and thermal loads on buckling of large-diameter cylinders typical of those 
which might be considered for an interstage structure of a launch vehicle and a fuselage 
of a supersonic transport. 


SYMBOLS 

The units used for the physical quantities in this report are given both in U.S. 
Customary Units and in the International System of Units (SI). A table of conversion fac- 
tors is given in appendix A. The relationship between these two systems of units can be 
found in reference 18. 

A cross-sectional area of stiffener 


Ai,Bi,Ci 4 X 4 matrices defined after equation (C16) 


'ij 


D 


1 ] 


E 


extensional stiffness of stiffened shell wall (eqs. (BIO)) 
bending stiffness of composite shell wall (eqs. (BIO)) 
Young’s modulus 


Ex,Ey Young's moduli of orthotropic layer of shell 


2 



G 


®xy 

I 


J 

Kij 


M 


MrpxjMxy 

Mxs^^Tr 

N 

Nx>Ny,Nxy 


Nx 

NxxjNxy 

Nxs,Nxr 


R 

T 


T 


cr 


u,v,w 


z 


shear modulus 

shear modulus of orthotropic layer of shell 
moment of inertia of stiffener about its centroid 
torsional constant for stiffener 

stiffnesses associated with coupling between bending and extension 
(eqs. (BIO)) 

buckling-moment variation with x as defined by equations (C 13) 


moment resultants for stiffened shell 


thermal moment resultants in multilayered shell 


thermal moment resultants in stringers and rings, respectively 


number of rings 


stress resultants for stiffened shell 


magnitude of applied compressive load at ends of cylinder 


thermal stress resultants in multilayered shell 


thermal stress resultants in stringers and rings, respectively 


radius of cylinder (to reference surface) 


change in temperature 


skin temperature change at buckling 


buckling-displacement functions of x defined by equations (C 13) 


curvature parameter for ring-stiffened cylinder bay, 



Ni= 


I^A 


Zi = 


u 

V 

1 w 

iM 


length of cylinder 


stringer spacing 


distance from reference surface of cylinder to surface on which Nx acts 


function 


shape factor for exponential temperature variation (eqs. (B5) to (B7)) 


height of stringer or ring 


integer 


total number of finite-difference intervals along length of cylinder 


ring spacing 


n 


number of full waves in cylinder buckling pattern in circumferential direction 


P 

Po 

t 

f 


pressure (positive for external pressure) 


sea-level atmospheric pressure 


thickness of isotropic cylindrical shell 


effective wall thickness of stiffened isotropic cylinder, — ^ + t 

d 


u,v,w displacements of stiffened shell at reference surface in x-, y-, and 

z -direction, respectively 
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x,y,z orthogonal curvilinear coordinates (see fig. 1) 

z distance from centroid of stiffener to reference surface (see fig. 1), positive 

if stiffener lies on external side of reference surface 

Aij structural coefficients defined after equations (Cl 5) 

Ilg strain energy of composite shell wall (shell plus stiffeners) 

IIl potential energy of loading 

a coefficient of linear thermal expansion 

o^,ay coefficients of linear thermal expansion for orthotropic layer of shell 

^sk>^r>^s temperature distribution factors (see eq. (Cl)) 

/3i,^2»* • ‘^7 structural coefficients in equations (C4) and (C6) 

^sk’^r>^s temperature distribution factors (see eq. (Cl)) 

shearing strain 

A distance between adjacent finite-difference stations 

poO 

6(x - ji) Dirac delta function defined so that \ f(x) 6(x - jl) dx s= f(jZ) where 

-OO 

6(x - ji) = 0 when x # jl 

®iH Kronecker delta, equal to zero when i H and equal to 1 when i = H 

ex,Cy normal strains 

changes in curvature 

jLt Poisson's ratio for isotropic material 

MxjMy Poisson's ratios for extension of orthotropic layer of shell 
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normal stress 


a 

T shear stress 

Subscripts: 

A prebuckling state 

B small changes away from prebuckling state which occur at buckling 

H . integer denoting finite-difference station at which ring is located 

cr at buckling 

i index of finite -difference station 

max maximum 

r rings (circumferential stiffening) 

s ‘stringers (longitudinal stiffening) 

sk skin 

Superscripts: 

i layer index 

t total strain in composite shell 

A subscript preceded by a comma indicates partial differentiation with respect to 
the subscript. 

Primes denote total derivatives with respect to x. 


6 



BUCKLING THEORY 


Assumptions 

The theory employed herein is an extension of the theory of reference 17. The 
cylinder considered (see fig. 1) is composed of a multilayered orthotropic shell, stiffened 
by uniform, equally spaced rings and stringers, all having linearly elastic properties. 

The elastic constants of the multilayered shell are taken as those given in reference 19. 
Individual layers of the shell are assumed to be such that their principal axes of orthot- 
ropy are alined with the principal directions of loading. The stringers are assumed to be 
closely spaced so that their elastic properties may be averaged over the stringer spacing. 
The effects of the rings are not averaged but instead the rings are considered to be located 
along the length of the shell at uniform spacings. The effects of the individual (discrete) 
rings are accounted for through the introduction of the Dirac delta or impulse function as 
described in references 17 and 20. The cylinder is assumed to be loaded by compressive 
and/or pressure loadings (mechanical loading), and is subjected to thermal stresses. 

In considering thermal stresses and deformations, the usual assumptions of engi- 
neering thermal-stress studies are adopted in that coupling between heat transfer and 
deformation is ignored. Thermal effects are accounted for by specifying axisymmetric 
temperature distributions in the shell and stiffening elements. The theory is sufficiently 
general to include arbitrary temperature variations through the thickness of the shell and 
depth of the stiffening elements. The effects of axisymmetric temperature variations 
along the length of the shell are accounted for in an approximate manner. 

The theory was derived by specifying strain energy expressions corresponding to 
nonlinear Donnell-von Karman displacements in the shell and stiffeners, and by applying 
the method of minimum potential energy to obtain equilibrium equations and associated 
boundary conditions. These equations were separated by a perturbation procedure into 
equations governing axisymmetric behavior prior to buckling and equations governing 
behavior at buckling. A complete development of the theory is contained in appendix B. 
Since similar developments have appeared in the literature (for example, refs. 17 and 21), 
only pertinent equations are presented. 

In the derivation of the governing equations, the ends of the stiffened cylinder were 
considered to be free to expand longitudinally (u 0). With this assumption, thermal 
buckling is a consequence of circumferential compressive stresses introduced by radial 
restraint at the boundaries (w = 0) and restraints resulting from differences in expansion 
between the stiffeners and the shell. 
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Governing Equations 


Prebuckling .- The equation governing the prebuckling behavior of a multilayered 
stiffened cylinder subject to axisymmetric thermal and mechanical loadings derived from 
the present theory (see appendix B (eq. (B21)) can be written as 


Dll 

\ 



2 1 ^ 11^12 
r \ Cii 







R 


^Ty ^Tr 


§^(Nx - Nts - Ntx) 


+ N-r(x) - p 


( 1 ) 


The equation is a Linear, fourth-order, ordinary differential equation in w^, the radial 
displacement of the cylinder prior to buckling. The coefficients of involve the 

structural stiffnesses of the composite shell wall (including stiffeners) Dij, C^j, and 
Kij which are associated with bending, extension, and coupling between extension and 
bending, respectively. Because rings are treated as discrete members, the circumfer- 
ential extensional stiffness C 22 is a function of x, the axial coordinate. Thus, equa- 
tion (1) contains involved variable coefficients and requires numerical techniques for 
solution. The magnitude of the applied axial compressive load Nx appears on the left- 
hand side of equation (1) as NxW^, a term which represents the contribution of nonlinear 
von Karman strain-displacement terms to the prebuckling analysis. 

The right-hand side of equation (1) involves loading terms resulting from tempera- 
ture gradients in the stiffened shell structure (terms with subscript T) and the applied 
external pressure p. The thermal loads Nipjjj Nrpy, Nipg, and are results of 

temperature gradients through the thickness of the shell and depth of the stiffening ele- 
ments, whereas N-p(x) contains terms associated with exponential variations of temper- 
ature along the length of the stiffened shell. The effect of thermal gradients and mechan- 
ical loadings is both to stress the cylinder and change its shape prior to buckling. The 
total effect can be accounted for by finding w^ and its derivatives by solving equa- 
tion (1) for suitable boundary conditions and applied loadings. 

Buckling . - The equations governing the buckling behavior of a stiffened shell 
derived from the present theory (see eqs. (B22)) can be written as 



(2a) 


(2b) 
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(2c) 


^ ^^B,xy ■ ^y^Bpcy ' ®^yB,yy 


+ ^xWjB^ - NxgW^ 


- %A^B,yy " 2NxyAWB,xy = « 

Equations (2a) and (2b) involve the stress resultants associated with extension and 


and shear 


during buckling. Equation (2c) involves the moment resultants 


^Mxg, Myg> ^yxgj ^xyg^ 3.nd contains some of the effects of deformations and 

stresses in the axisymmetric prebuckling state (terms with subscript A). In form, 

equations (2) reduce to more familiar classical buckling equations if the term Nx wV 

B " 

is omitted in equation (2c). The classical buckling theory contains the effects of pre- 
buckling loads ^Nx, Ny^, and Nxy^j but not prebuckling deformations ^w^ and 

The effects of thermal gradients appear in equations (2) in two ways. First, the 
buckling stress and moment resultants are functions of w^ so that terms corresponding 
to thermal prebuckling deformations appear in equations (2a), (2b), and (2c). Secondly, a 
direct thermal-stress term appears as Ny^Wg in equation (2c). This term repre- 
sents the effects of circumferential stresses introduced by differences in expansion due to 
thermal gradients or variation in thermal properties in the stiffened shell wall. In many 
existing theories for thermal buckling, the direct thermal-stress term is retained, 
whereas the less obvious effects of thermal deformations are omitted. 


Numerical Solutions 

The complexity of equations (1) and (2) requires the use of numerical techniques. 
Solutions to the equations governing prebuckling as well as buckling behavior were 
obtained by omitting the applied torsional loadings ^^xy^^ and employing finite differences. 

Both sets of equations were formulated as a system of second-order difference equations 
which are documented in appendix C. The prebuckling equations (eqs. (B18) and (B21)) 
were solved by matrix algebra by exploiting Gaussian elimination. Buckling loads were 
extracted from the equations governing buckling (eqs. (CIS) and (C16)) by either of two 
methods: determinant plotting or modal iteration. In the determinant plotting solution, 
either the applied axial compressive load N^x> the pressure p or a base temperature 
T^ can be used as a buckling parameter. For the modal iteration solution, T^ must 
be the buckling parameter. The parameter T^ is related to temperature changes 
between the cylinder shell and stiffeners and the boundaries of the stiffened cylinder. 
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The two numerical solutions have been incorporated into computer programs 
entitled BAMSOC I and BAMSOC II (Buckling Analysis of Multilayered Stiffened 
Orthotropic Cylinders), The programs were written in FORTRAN IV for the Control 
Data 6600 computer. Input variables and sample problems are given in appendix D. The 
cards and computer listings for these programs can be obtained from COSMIC, University 
of Georgia, Athens, Georgia 30601. 

APPLICATIONS OF ANALYSIS 


By use of the present analysis, the thermal buckling behavior of four types of 
cylindrical shells was investigated: unstiffened, ring stiffened, stringer stiffened, and 
ring and stringer stiffened. Unstiffened and ring-stiffened shells were studied in order 
to relate results from the present theory to published results from less complex theories. 
The behavior of two large-diameter shells — one stringer stiffened, the other ring and 
stringer stiffened - was also studied. The interaction of thermal and mechanical loading 
in these two structures was investigated in order to assess the sensitivity of contempo- 
rary aerospace designs to thermal buckling. 


Unstiffened Cylinders 

A solution for buckling of unstiffened cylinders subjected to uniform temperature 
has been presented by Hoff in reference 1. In the reference, calculations of the thermal 
buckling temperature were made for a steel cylinder with a radius-thickness ratio of 300 
and a length-radius ratio of about 0.3. The boundary conditions employed corresponded 
to those of simple support = Vg = Wg = = oj . The cylinder was found to buckle 

at an unusually high temperature. 

In calculations using the present theory, both the determinant plotting solution and 
the modal iteration solution were employed in an attempt to find a buckling temperature 
for the cylinder. The modal iteration solution would not converge. Further investiga- 
tion was conducted with the determinant plotting solution, but no buckling temperature 
could be found within a practical temperature range. 


The basic difference between the present theory and that of reference 1 lies in the 
treatment of the effects of deformations and stresses in the shell prior to buckling. In 
the theory of reference 1, the effects of thermal stresses introduced by restraints offered 
by the boundaries of the cylinder are accounted for by introducing circumferential thermal 
stresses into Donnell's buckling equation. In the Donnell equation, however, the effects of 
rotations and changes in rotations ^w^ and w^^ in the cylinder prior to buckling are 


ignored. The present theory accounts for thermal effects in a more consistent fashion, 
and as previously discussed, contains additional terms in the buckling equations which 
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are associated with shell prebuckling rotations. In order to make computations for the 
steel cylinder on a more comparable basis, these additional terms were suppressed in 
buckling calculations. Thermal buckling predictions were then in excellent agreement 
with the results of reference 1. It must be concluded that the additional terms of the 
present theory tend to alleviate the effects of the boundary restraints; thus, the buckling 
temperature of this cylinder is beyond a practical range of interest. 


An approximate solution for buckling of clamped, unstiffened cylinders subjected to 
uniform temperature is presented by Sunakawa in reference 6. The solution is based on 
a theory which includes the prebuckling displacement effect just discussed. The clamped 
boundary conditions for buckling implied in the solution are Njjg = 0, Vg - 0, Wg - 0, 


and 


Wg^ = 0. Thermal buckling predictions for the steel cylinder previously described 


were performed by using the same clamped boundary conditions in the present theory. 
Results for the nondimensional buckling temperature difference oTcr and the number of 
circumferential waves n in the buckling mode are as follows: 



Present 

theory 

Theory of 
reference 6 

oTcr 

n 

0.01440 

46 

0.01255 

44 


It can be seen that the present theory yields a buckling temperature which is about 15 per- 
cent higher than that predicted by reference 6 with similar circumferential mode shapes. 
The differences in buckling predictions are attributed to lack of convergence of the one- 
term Galerkin solution of reference 6. Convergence of the present solution was investi- 
gated by varying the number of finite-difference stations and was found to be satisfactory 
when 100 finite-difference stations were employed. 

Ring-Stiffened Cylinders 

A theoretical model for buckling of ring-stiffened cylinders has been proposed by 
Anderson in reference 10. In Anderson’s theory, the ring-stiffened cylinder is idealized 
by considering the behavior of a skin bay in the central regions of an infinitely long cyl- 
inder. Thermal stresses caused by temperature differences between the skin and rings 
are found for the bay by assuming that no rotation occurs at the rings. The circumfer- 
ential stresses are then introduced into Donnell's buckling equation for an unstiffened, 
isotropic shell. Boundary conditions (at the rings) during buckling are assumed to be 
those corresponding to conventional simply supported or clamped boundaries. 

In contrast, the present theory treats a cylinder of finite length. Thermal stresses 
and deformations prior to buckling are found for a cylinder supported at its ends in 
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conventional fashion, but with intermediate, discrete stiffening rii^s. Rotation of the 
rings may occur during thermal loading. Thermal stresses and deformations for a dis- 
cretely stiffened shell are introduced consistently into the buckling equation. Boundary 
conditions during buckling are enforced only at the ends of the cylinder so that the con- 
tinuity conditions at each ring are natural consequences of the structural stiffnesses of 
the ring and shell at the line of attachment. 


Comparison with reference 10 .- In order to compare Anderson’s predictions with 
the present results, an analytical investigation was made of the thermal buckling behavior 
of ring-stiffened cylinders with proportions similar to the test cylinders of reference 11. 
These cylinders were internally stiffened with equally spaced Z -shaped rings having the 
dimensions shown in figure 2, Geometric and mechanical properties employed in the 
present calculations are indicated in table I, In the investigation, the rings and bound- 
aries of the cylinder were considered to be at room temperature, whereas the cylinder 
skin was heated to a constant temperature. The nonuniform temperature in the stiffened 
cylinder produces circumferential thermal stresses that cause buckling of the cylinder 
skin between rings. The cylinder was assumed to be simply supported at its ends with 
the boundary conditions Nxg = 0, Vg = 0, Wg = 0, and Mxg = 0. 


For the comparison, the length, radius, and thickness of the cylinder were held con- 
stant for all the computations. The number of rings in the cylinder (and hence the ring 


spacing 1 ) was varied to study the effects of the bay curvature parameter 



on thermal buckling characteristics. The results of reference 10 suggest that somewhat 
unusual behavior occurs for small Z in ring-stiffened cylinders. The overall cylinder 
length a was much shorter than that of the test cylinders of reference 11. The particu- 
lar value selected (see fig. 2) was based on considerations of the number of rings and 
finite-difference stations required to investigate small values of Z. 


Results of the theoretical buckling predictions from reference 10 and the present 
theory are compared in figure 3. The nondimensional skin-buckling temperature oTcr 
is shown as a function of Z. Results from the present theory are shown as discrete 
points corresponding to calculations with a finite number of rings. The solid curve has 
been faired through the points to indicate trends. The dashed curves on the figure are 
based on results presented by Anderson (ref. 10) for simply supported and clamped cylin- 
der bays with rigid rings, A correction which accounts for flexibility of the rings is sug- 
gested in the reference. When the rings are considered to be flexible, circumferential 
thermal stresses near the rings are less than those computed for rigid rings; therefore, 
the thermal buckling temperatures are higher. For the present cylinders the correction 
amounts to an increase of about 7 to 7^ percent above the values shown in figure 3. 
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In figure 3, it can be seen that results from the present theory are in reasonable 
agreement with Anderson’s results for clamped cylinder bays when Z is greater than 
15. The numerical differences shown for this range of Z are about 5 percent. For 
smaller Z, a larger disparity exists. Studies of prebuckling stress distributions based 
on the present theory suggest that the idealizations used by Anderson to obtain thermal 
stresses are reasonably accurate. It is believed, therefore, that the analytical differ- 
ences shown in figure 3 stem mainly from the additional thermal deformation terms in 
the present buckling theory, and the effective ring boundary conditions during buckling. 

The present results suggest that when Z is greater than 15, the cylinder rings offer 
nearly clamped support to the skin bays during buckling. Furthermore, both the present 
theory and that of reference 10 suggest that the thermal buckling temperature becomes 
independent of ring spacing for this range of Z. The theoretical differences and unusual 
behavior for cylinders with smaller Z is probably associated with the interaction of 
closely spaced rings. 

Thermal-mechanical interactions .- A further topic of interest in the thermal buck- 
ling behavior of ring-stiffened cylinders is the interaction between thermal and mechan- 
ical loading. The behavior of the simply supported cylinders just discussed was investi- 
gated under combinations of axial compression and thermal loads induced by heating the 
cylinder skin. Results of interaction studies are presented in figure 4, where the nondi- 
mensional mechanical buckling load has been plotted against the thermal buck- 

ling load oTcr* Buckling interaction curves are presented for cylinders of fixed overall 
geometry (see fig. 2) with three, five, and seven rings. The buckling loads shown for 
axial compression loading alone ^oTcr = o) were compared with published results for iso- 
tropic cylinders having a similar value of Z (ref. 22). The comparison suggested that 
the rings of the present cylinders provide almost clamped supports for the skin bays 
when loaded in axial compression. A similar result has already been noted for thermal 
loading when Z is large. The interaction curves shown are very regular in shape with 
the exception of the case of three rings. Results of the interaction curves were compared 
with interactions suggested by reference 10 for clamped cylinders but no correlation was 
evident for the range of curvature parameter considered (Z = 8 to 30). In this range of 
Z, the results of reference 10 suggest that substantially larger thermal loads can be car- 
ried when axial compressive loads are present, whereas results from the present theory 
do not support this conclusion. 


Stringer-Stiffened Cylinder 

In order to investigate the thermal buckling behavior of longitudinally (stringer-) 
stiffened shells, the behavior of a large-diameter aluminum shell representative of pro- 
portions which might be considered for an interstage structure in a large launch vehicle 
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was investigated. The structure considered is an integrally stiffened shell which has 
been used in theoretical and experimental studies of stiffener and loading eccentricity 
effects. (See refs. 23 and 24.) The geometry and mechanical properties adopted for the 
present investigation are given in figure 5 and table 1(b). Thermal buckling problems in 
interstage structures might conceivably arise from thermal gradients induced by rapid 
aerodynamic heating of the shell outer skin; as a result, thermal lags between stiffener 
and shell temperatures occur. For the present analysis, the skin was assumed to be 
heated whereas the stringers remained unheated. Bending and thrust loads are also 
present during flight; thus, it seemed desirable to investigate the interaction between 
axial compressive and thermal loads. For purposes of the investigation, the cylinder 
was assumed to be clamped at its ends ^Nxg = Vg = Wg = Wg ^ = oj. Local buckling of 

the skin was not considered. 

Prebuckling stress d istribution.- The stress distributions induced by combinations 
of axial compression and temperature were investigated by use of the results of the solu- 
tion of the equations governing deformations and stresses prior to buckling. (See eq. (1).) 
Stress distributions in the skin and stringers of the interstage for three combinations of 
applied compressive load Nx and change in skin temperature Tg^ are presented in 
figure 6 as functions of the distance x/a measured from the supported end of the shell. 
In figure 6(a) the axial stress at the skin middle surface is shown to increase in propor- 
tion to the applied axial load. The axial stress asymptotically approaches the mem- 
brane state of stress = Njjy/Et j in the central areas of the shell. The higher values of 
stress near the ends of the shell (x/a = 0) are due to bending stresses developed from end 
constraints. In figure 6(b), the circumferential stresses at the skin middle surface shown 
increase roughly in proportion to skin-temperature changes. The circumferential 
stresses are small in the central regions of the cylinder and the stresses shown near the 
ends are largely due to restraints in expansion between the cold boundaries and cold 
stringers and the heated shell. A significant difference between circumferential stress 
variation in the present stiffened shell and results for unstiffened shells of similar geom- 
etry is that stresses decay much more rapidly in unstiffened shells. (See ref. 13.) 

Stresses developed in the longitudinal stiffening for various combinations of axial 
compression and thermal load are shown in figure 6(c). The stress distribution in the 
stringer is plotted against the normal coordinate z nondimensionalized by the stringer 
height h for various longitudinal stations x/a. The origin of z was taken to be the 
surface of the stiffener in contact with the shell. The stress distributions shown in fig- 
ure 6(c) indicate that bending stresses in the shell increase appreciably as the tempera- 
ture difference between skin and stringers increases. For the case with no applied axial 
compression, the calculated stress in the outstanding flange of the stiffener actually 
exceeds the ultimate tensile strength of aluminum. 


14 



The results presented for the stress distribution in the interstage structure suggest 
that complex states of stress exist which depend on the relative magnitude of the applied 
axial compression and the differences in expansion between the boundaries and stringers, 
and the heated shell. Studies of the detailed stress distribution in skin and stiffeners 
indicate that intolerably large bending stresses exist in the stringers near the ends when 
only thermal loads are present. The results also indicate that thermal stresses are 
developed from boundary effects which appear to persist over relatively large sections 
of the longitudinally stiffened shell. 

Buckling of stringer -stiffened cylinder .- The buckling behavior of the idealized 
interstage structure was investigated under combinations of axial compressive and ther- 
mal loads. Results of the investigation are presented in figure 7 as an interaction curve 
between the normalized axial compressive load at buckling Nx^p^Et and the change in 

skin temperature at buckling oTcj.. The three combinations of axial load and tempera- 
ture change discussed in the section entitled "Prebuckling stress distribution" are shown 
as cross marks in figure 7. The buckling load shown for the clamped shell under axial 
compression alone (oTcr ® o) was in good agreement with published results, the present 
calculation being about 2 percent higher than the result presented in reference 24. 

In figure 7 it can be seen that the interaction curve has an unusual shape when small 
thermal loads are present. For temperature changes of about 50° F (28 K) in the skin, 
that is, oTcr = 0.0007, the shell buckling load actually increases to about 1.4 times the 
compressive buckling load of an unheated shell. It appears that small thermal loads can 
introduce favorable deformations prior to buckling which increase the shell's effective 
stiffness to resist axial compression. As can be seen from the interaction curve, how- 
ever, when thermal loads became large, the longitudinally stiffened shell had little 
resistance to the circumferential thermal stresses so that buckling can occur by thermal 

= (^. The temperature change corresponding to this case is about 

130° F (72 K), and is well within a practical temperature range for aluminum. In inter- 
stage structures, when thermal environments are present, simple insulators such as cork 
are used to shield the outer skin of the structure. However, in many instances, the insu- 
lation is not adequate to eliminate all thermal gradients. The present results suggest 
that small gradients may be tolerated, but if large thermal gradients are present, the 
structure should be carefully considered with regard to the possibility of thermal buckling. 

The normal component of the buckling mode shape for the interstage structure under 
thermal loading alone is shown in figure 8, where the normalized variation of the radial 
buckling displacements w/Wmax plotted against the length parameter x/a where 
a is the length of the shell. For comparison, prebuckling deformations at the onset of 
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buckling are shown in the inset in figure 8. The buckling wave pattern is characterized 
by a nearly sinusoidal wave along the length of the shell and many sinusoidal waves 
around the circumference (n = 65). Because the number of circumferential waves is 
large, the nodes in the buckle pattern in the circumferential direction are separated by 
distances which are about 1-| times the stringer spacing. Since the stiffness properties 
of longitudinal stiffening are averaged in the present theory, the calculated circumferen- 
tial mode shape is one which is near the limits of applicability of the theory. 

Ring- and Stringer-Stiffened Cylinder 

In order to investigate the thermal buckling behavior of ring- and stringer-stiffened 
shells, the behavior of a stiffened cylinder with geometry similar to that which might be 
considered in a preliminary design of a supersonic transport was investigated. A rather 
long, uniform, titanium fuselage section, bounded forward by a heavy bulkhead and aft by 
the wing- joint structure, was selected. The geometry and mechanical properties adopted 
for the present calculations are presented in figure 9 and table I. Studies of design 
loads and predicted thermal environment for this type of structure suggested the possi- 
bility of thermal buckling as a consequence of a rapid climb to altitude. In this maneuver 
there is a thermal lag between the skin and stiffening elements so that for a short period 
of time the skin is at elevated temperatures while both rings and stringers are unheated. 
Bending and pressurization loads on the fuselage are also present during the maneuver; 
therefore, it seemed desirable to investigate the interaction of axial compressive loads 
and thermal loads for various internal pressures. For purposes of the investigation, the 
cylinder was assumed to be simply supported at its ends ^Nxg = Vg = Wg = = 0^ and 

the axial compression was introduced by means of a line load applied at the skin middle 
surface. 

Prebuckling stress distribution . - The stress distribution induced by combinations 
of axial compression and temperature was investigated by using results based on the 
solution of equation (1). Stress distributions in the skin, stringers, and rings of the fuse- 
lage under an applied compressive load of 1650 Ibf/in. (0.289 MN/m) and temperature 
difference of 360*^ F (200 K) between skin and stiffeners were obtained. Results are pre- 
sented in figure 10. The diffusion of load into each element of the structure is illustrated 
by showing results for various stations near one of the supported ends of the cylinder. 

In figure 10(a) the stress distribution in the sMn of the cylinder is shown as a function of 
length. At the end of the cylinder (x/a = 0), both axial and circumferential stresses are 
very high, but are still within the elastic range of titanium. The stresses are not nearly 
as large as those for the clamped interstage structure already discussed. Sharp peaks in 
stress occur at each of the 19 rings in the shell. The circumferential stress shown in 
bays that are distant from the boundary is a result of the differences in thermal expansion 


16 



I 


between the rings and the shell. The curves shown in figure 10(a) indicate that the state 
of stress does not vary from bay to bay beyond the third bay (x/a *=0.15). 

Variations in the stress distribution in the stiffening elements are shown in fig- 
ure 10(b). The assumptions of the present theory are such that the ring is in a state of 
uniform, uniaxial stress (in this case, tension). The ring stress variations shown in fig- 
ure 10(b) indicate that the interior ring stresses are essentially constant in rings that are 
not adjacent to the boundary. The assumptions of the present theory are such that the 
stringers (longitudinal stiffening) are also in a state of uniaxial stress, although bending 
stresses through the depth of the stringers may exist. The magnitudes of these bending 
stresses are given in figure 10(b) by plotting stresses against the depth parameter z/h 
for various longitudinal stations. The bending stresses in the interior of the shell are 
mainly due to thermal stresses and restraints to expansion offered by the rings. The 
bending stresses near the ends of the cylinder also include the effects of introducing the 
applied axial compression load by means of a line load on the skin midplane. The edge 
moment introduced by the eccentricity of the line load with respect to the stringer-shell 
centroid has been noted in previous studies (ref. 23) and large effects on buckling predic- 
tions have been found. In the present case, the stress distributions shown in figure 10(b) 
suggest that this edge moment decays rapidly from the cylinder ends so that the state of 
stress in the stringers does not vary from bay to bay beyond the second ring. 

The results presented for the stress distribution in this example indicate the com- 
plexity of the state of stress existing in a rii^- and stringer-stiffened shell under com- 
bined mechanical and thermal loadings. Studies of the detailed stress distributions in 
skin and stiffening elements suggest that no large stresses exist that are not within the 
capabilities of the structural material. The effects of the boundaries of the stiffened 
shell on stress distribution appear to be negligible beyond the second bay of the cylinder. 

Bucklin g of ring- and stringer-sti ffened cylinder .- The buckling behavior of the 
idealized supersonic transport design was investigated under combinations of axial com- 
pressive and thermal loads. In addition, the effects of including cabin pressurization as 
an applied load were also included. The results of the investigation are summarized in 
figure 11 by means of interaction curves between axial compressive load at buckling 
Nxcr/^^ and relative skin temperature at buckling cfUcv solid curves shown cor- 

respond to various values of internal pressure p normalized by sea-level atmospheric 
pressure Pq. The maximum value of p/pq shown is about that proposed for cruising 
altitude. 

Results shown in figure 11 for p/pq = 0 indicate that the interaction between axial 
compressive and thermal loads is substantial. In fact, in the absence of axial compres- 
sion, thermal loads are large enough to cause buckling in a practical temperature range. 
The value of cfTcr shown for this case corresponds to a temperature difference of about 
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525° F (292 K) between the skin and stiffening elements. Results shown in figure 11 for 

cases where internal pressurization exists (p/Pq < o) indicate that even small values of 

internal pressure are very beneficial in alleviating the problem of buckling under thermal 

loading alone. At cabin pressure at cruise (p/Pq = -0.82), the interaction curve is nearly 

horizontal over a practical temperature range and the buckling load is roughly 1— times 

2 

that of an unpressurized fuselage. In some aircraft design criteria for buckling, how- 
ever, the beneficial effects of pressure are not included in buckling calculations; there- 
fore, the fuselage structure is required to be stable even if all cabin pressure is acciden- 
tally lost. Under these circumstances, elimination of buckling under thermal and com- 
pressive loads or thermal loads alone may be a serious design consideration. 

An additional effect shown in figure 11 by the dashed curve is that of introducing the 
axial compressive load at the neutral surface of the sheet-stringer combination 
(e = 0.13 in. (0.33 cm)) rather than at the skin midplane (e = 0). The dashed curve sug- 
gests substantial improvement in buckling strength under axial compression alone due to 
the edge moment effect already mentioned. However, the interaction curve is unaffected 
when large thermal loadings are present. 

The normal component of the buckling mode shape for the stiffened shell under 
thermal loading alone with no internal pressure is presented in figure 12, where the nor- 
malized variation of the radial buckling displacements plotted against 

length x/a for one-half of the cylinder. For comparison, prebuckling deformations at 
the onset of buckling~are shown in the inset in figure 12. The wave shape shown in fig- 
ure 12 is typical of a panel instability buckling pattern in that nodes occur at each of the 
rings. It can be seen that the amplitude of W decays in the interior bays so that in the 
central bays of the cylinder almost negligible deformation is present. The variation in 
amplitude from bay to bay suggests that buckling is a consequence of restraints at or 
near the boundaries of the cylinder. The presence of longitudinal stiffening apparently 
makes thermal edge effects persist over a large length as evidenced by the large ampli- 
tude of the buckle wave in the fifth bay of the cylinder (0.20 ~x/a ~ 0.25). As was the 
case with the interstage structure, a large number of circumferential waves (n = 57) 
occurred in the buckling mode when thermal loading predominated. In this case, however, 
at least three stringers are included in each circumferential half-wave. Thus, the buck- 
ling mode found appears to be a legitimate panel instability mode and is consistent with 
theory employing the average stiffness properties of longitudinal stiffening. 

The present study of an idealized supersonic transport design indicates that a poten- 
tial design problem exists in that temperature gradients which might be induced by a 
rapid climb to altitude were sufficient to buckle the shell in a practical temperature range. 
The presence of internal cabin pressure is sufficient to alleviate the thermal loads. If 
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internal cabin pressure is not available as a stabilizing force, additional reinforcement 
must be judiciously applied to the structure. 

CONCLUDING REMARKS 

An extension of the theory of NASA Technical Note D-4283 for buckling of stiffened 
cylindrical shells under combinations of mechanical and thermal loadings has been pre- 
sented. It includes the effects of eccentricity in multilayered reinforcement, deforma- 
tions prior to buckling, and discrete rings. The theory is sufficiently general to account 
for the averaged stiffnesses of stringers and multiple orthotropic shell layers as well as 
arbitrary temperature distributions through the thickness of the shell and depth of the 
stiffeners. Buckling solutions have been obtained by employing finite differences. Two 
computer programs from which buckling loads can be found (either by determinant plotting 
or modal iteration) have been developed. 

Theoretical results for thermal buckling have been compared with existing solutions 
for unstiffened and ring-stiffened cylindrical shells. Results from the present theory for 
a simply supported unstiffened shell indicated that thermal buckling would not occur in a 
practical temperature range. Results for a clamped unstiffened shell were in agreement 
with published theoretical results. Buckling predictions for a ring-stiffened shell with 
various ring spacings were in agreement with existing theoretical results when the cylin- 
der bay curvature parameter is large. In cases where large disparities between pub- 
lished theoretical results and those presented herein exist, differences are attributed to 
the effects of prebuckling deformations which have not been accounted for in previous 
theories. 

The effects of interaction of thermal loading and axial compression on the buckling 
behavior of two types of contemporary stiffened shell structures have been investigated. 
Thermal loads were introduced by considering the skin of the structure to be heated, 
while the stiffening elements and shell boundaries remain unheated. The behavior of an 
aluminum, clamped, large-diameter longitudinally stiffened cylinder representative of a 
launch vehicle interstage structure was investigated. Stress distribution studies indi- 
cated that thermal stresses induced by moderate skin temperatures cause large bending 
stresses in the stiffeners near the ends. Buckling calculations suggest that the inter- 
stage structure may buckle under thermal loading alone at a realistic skin temperature. 

The behavior of a titanium simply supported, ring- and stringer-stiffened cylinder 
representative of a preliminary design for a supersonic transport fuselage was also 
investigated. Stress distribution studies indicated that no large thermal stresses beyond 
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the capabilities of the shell material were present. Buckling calculations suggest that 
the fuselage also can buckle under thermal loading alone at a practical skin temperature. 

Near the completion of this work, preliminary calculations for the ring-stiffened 
cylinders by another investigator suggested that thermal-buckling behavior could be 
influenced by the type of ring theory employed in the analysis. For the ring- stiffened 
cylinders, out-of-plane bending stiffness appeared to be an important parameter. To 
account for this effect, a ring theory which includes out-of-plane bending behavior was 
developed and is presented in appendix E. This ring theory was incorporated into one of 
the computer programs; however, no numerical results were available for this paper. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va. , April 28, 1970. 
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APPENDIX A 


CONVERSION OF U.S. CUSTOMARY UNITS TO SI UNITS 

The International System of Units (SI) was adopted by the Eleventh General Con- 
ference on Weights and Measures in 1960. (See ref. 18.) Conversion factors for the 
units used in this report are given in the following table: 


Physical quantity 

U.S. Customary 
Unit 

Conversion factor 
(♦) 

SI Unit 
(**) 

Length 

in. 

0.0254 

meters (m) 

Stress modulus 

ksi 

6.895 X 106 

newtons/meter2 (N/m2) 

Stress resultant 

Ibf/in. 

175.1 

newtons/meter (N/m) 

Temperature change .... 

1 _ « 

Op 

5/9 

Kelvin (K) 


^Multiply value given in U.S. Customary Unit by conversion factor to obtain 


equivalent value in SI Units. 

**Prefixes to indicate multiple of units are as follows: 


Prefix 

Multiple 

giga 

(G) 

109 

mega 

(M) 

106 

kilo 

(k) 

103 

deci 

(d) 

10-1 

cent! 

(c) 

10-2 

mi Hi 

(m) 

O 

1 

GO 
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APPENDIX B 


DEVELOPMENT OF GOVERNING EQUATIONS 


Stress-Strain Relations 


In the ith layer of the shell, the stress-strain relations can be written as 




K 


ri[4 

n * ^ 


My4 


MxMy 


E. 


1 ~ UL^ u} 






= G^ 
xy^xy 




(Bl) 


where and oy are coefficients of linear expansion associated with thermal elonga- 
tion in the axial and circumferential directions, respectively, T is a change in temper- 
ature relative to a thermal-stress -free state and the remaining quantities with the super- 
script i are orthotropic elastic constants associated with extension. The strains 
appearing in equation (Bl) are related to the shell displacements u, v, and w by 
Donnell-von Karman relations having the form 


with 




^y ^ ^ 

7xy - 7xy ^^xy 


(B2) 


^x =Upc 



2 






2 

,y 


>'xy = ",y + ^,x + ^,x^,y 


% = -w,xx 
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APPENDIX B — Continued 


'^y “ "^,yy 

/C^y = -2w xy 


In the stringers 



whereas in the rings 

<^Yr = Er^4 ' 

To account for axisymmetric temperature variation along the length of the shell, the 
relative temperature of the shell is expressed as 


T = Tsk 1 + Vsk®’ 


gx/a\ 


(B5) 


where a is the length of the shell, and and g are constants. Similarly, in the 

stringers. 


T = Tc 


+ roe 


gx/a 


(B6) 


and in the rings 



+ 



(B7) 


where Tgj^, Tg, and Tj. are functions of z. By specifying temperature distributions 
with equations (B5) to (B7), the first-order effects of longitudinal variation in thermal 
distributions can be obtained. The present theory, however, does not consider the effects 
of variation of mechanical properties with the assumed longitudinal temperature 
distribution. 


Potential Energy 

The integration of the familiar strain energy expression 




dV 
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APPENDIX B - Continued 


over the depth of the composite shell wall (shell plus stiffeners) yields the following 
expression: 

Hs = I < Nx€x + Ny6y + NxyVxy + MxKx + MyKy + (Mxy - Myx)-^ - |(Ntx + Nts)^ + (^Tx + Mxs)«x 

Surface L 

+ ^Nxy + Nxr)«y + (^xy + Mxr^ dx dy (BS; 


where 




~ ^11% ^12 ^y " ^^Tx ^Ts) 

Ny = Ci2^x + ^22^y *^12% *^22^y “ ^^Ty ^Tr^ 

Nxy = CQQ-y^ + Kg0Kxy 

Mx = + Ki26y + Dn/c^ + Dl2«y " (^tx + 

My = Kj2ex + K22Cy + ^12% + ^22>^y ~ (^Ty + ’^Tr) 

/ GgJg \ 

Mxy = "^ee'J'xy ^66%y + “2;^ '^xy/ 

N 

Myx = ^eerxy + ^^66%y + I Z GrJr%y 


j=l 


(B9) 




where 6(x - jf) is the Dirac delta function associated with the N rings on the cylinder. 

The quantities Cy, Djj, and Kjj appearing in equations (B9) are structural stiff- 
nesses associated with extension, bending, and extension-bending coupling, respectively, 
and are given by the following integrals: 


Cii = 


C12 = 


C22 = 


1 - 






^X^y 

1 - PxMy 

1 - 


dz 


P-xl^y 


dz + 


N 

^ 6(x - ji) 


ErA 


r"r 


j=l 


'66 




Gxy dz 


(BlOa) 


(BlOb) 


(BlOc) 


(BlOd) 
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*^11 

Ki2 

K22 

K66 

Dll 


1 - 


z dz + 


PxMy 


j 

— i — z dz 


1 - 


MxMy 


N 


£j 

^ z dz + / 6(x - iZ) ErAyZ^. 

j 1 - l^xi^y jfi 

y GxyZ dz 

^2dz + 


1 — jj.^jXy 


Di2 


z2dz 

1 - MxMy 



(BlOe) 

(BlOf) 

(BlOg) 

(BlOh) 

(BlOi) 

(BlOj) 

(BlOk) 

(BlOl) 


where the integration is to be performed over the multiple layers of the shell wall. A 
convenient numerical scheme for the integration for a specified reference surface can be 
found in equations (83) and (85) of reference 19 (p. 34). The properties associated with 
the stiffeners (location of the centroid of the stiffener z, the moment of inertia I, the 
cross-sectional area A, and the torsional constant J) are also defined (when appropri- 
ate) with respect to an arbitrary reference surface. In specifying the structural stiff- 
nesses, the dependence of the elastic constants upon temperature can be accounted for in 
the shell by analytically dividing the multilayered shell into thin layers and performing 
the required integration of equations (BIO) numerically. In the stiffeners, the following 
sample definitions can be employed to perform a comparable integration through the depth 
of the stiffeners: 
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= Es(z)dAs 

d d JAs 

=1 r Es(z)zdAs 
d d JAs 


The thermal moment and thermal force terms appearing in equa- 

tions (B9) are 


Ntx = 

Nxy = C^y^l + 

Nts = ^Ts(^ + 

N . 

^Tr ~ ^ ^Trr ■*■ y 

j=l ^ J 

Mtx = Ktx(^ ^sk®^'^^) 

Mxy = Kry(l + ygk®^^^) 

Mts = Kts(i + ^s®^^') 

M^r = ^ 6(x - jZ) Kpr(l + 


(Bll) 


(B12) 
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where 


ctx - 1 + ''y‘v)'rrj^ ’’sk 

Crps = ^ y “s^sTs 
As 


^Tr — \ 

A„ 


dAr 


Kj-x 

K-Ty 




Vy 


I (°y ^ Mitt 


Mx^y 


Tsk"= 


Kps “s^sTs^i dAg 

a * Ag 

Xrp^ “ ^ o^EpTpZ dAj* 


For Cipjj, C-py, Kjpx> and Kpy, the integration is to be performed over the 
multiple layers of the shell wall and can be accomplished by the method suggested for 
equations (BIO). 

As in reference 17 the potential energy of mechanical loading for the case of an 
applied compression load Nx, an applied shear load Nxy, and an external pressure p 
is given by 

p27rR c'B. /^27tR ^ . \i a p2ttR .. la 

nL = ^„ Jo P* '^y + Jo Mo^’' 


where e is the distance from the arbitrary reference surface to the line on which the 
applied mechanical load resultant Nx acts. 
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Equilibrium Equations and Boundary Conditions 


By the principle of minimum potential energy, the vanishing of the first variation of 
the total potential energy + 11^^ = o) yields the equilibrium equations and appropriate 

boundary conditions. The first variation with respect to the displacement variables u, 

V, and w yields 


Nx,x + ^xy,y - 0 


Nxy,x + Ny,y - 0 


) (B14) 


/ \ Ny 

+ (Mjjy - MyxjjXy " ^y,yy " NxW^xx " %^,yy " ^NxyW^xy + P = 0 


with the following boundary conditions to be specified at each end of the cylinder: 


Nx + Nx = 0 


Nxy + Njjy - 0 




Mx + Nx6 = 0 

^x,x “ ^^xy,y “ ^yx,yj ^x^^x ^xyW^y = 0 


or 


u = 0 
V = 0 

w,x = 0 
w = 0 


(B15) 


The nonlinear equilibrium equations and associated boundary conditions (eqs. (B14) 
and (B15)) are used to obtain equations governing the prebuckling and buckling states of 
the stiffened shell by taking the displacements u, v, and w to be 

u(x,y) = u^(x) + UB(x,y) 

v(x,y) = v^(x) + Vg(x,y) > (B16) 

w(x,y) = w^(x) + Wg(x,y) 

J 

In equations (B16) the subscript A denotes the axisymmetric, prebuckling displacements 
of the stiffened shell, and the subscript B denotes the infinitesimal nonaxisymmetric 
displacements that occur at buckling. 


Prebuckling Equations 

For axisymmetric prebuckling deformations, the equilibrium equations are found 
from equations (B14) by discarding terms which have derivatives of y. The resulting 
equations are 
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- » 

(B17a) 



(BlTb) 


k 

> 

> 

1 

> 

+ 

II 

O 

(B17c) 

with boundary conditions 

Nxa + Nx = 0 or u^ = 0 

(B18a) 


Nxy^ + Nxy = 0 or = 0 

(B18b) 


^x^ + Nxe = 0 or w^^^ = ° 

(B18c) 


^^Xa,x ^^a"^A,x = 0 ^A = 0 

(B18d) 

where the subscript A terms are defined by inserting u = Uj^(x), v = v^(x), and 
w = w^(x) in the definitions of the stress resultants (eqs. (B2) and (B9)). 


As discussed in the text, the cylinder was considered to be free to expand longitudi- 
nally ^u^ o) so that integration of equation (B17a) yields 


II 

(B19) 

where Nx is the applied compressive load at the boundary. From equation (B19) and 

the definition of Nv , 

^A 


Kii C,2W^ 

Cii c„ R di 

(B20) 

The use of the definitions of and Ny together with equation (B20) in equa 

A 

tion (B17c) yields the following equation in w^: 

- 



(B21) 
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with 


Nrrfx^ - + Crr 

^tW - ^Ftx + ^Ts ^2 ® ' 


- Kp: 


^sk^ ^gx/a ^ ^sg ^gx/a 


^ a2 


+ Kj. 


® a2 


where the primes denote total differentiation with respect to x. 


Buckling Equations 

The equilibrium equations and boundary equations which govern the buckling behav- 
ior of the stiffened cylinder are obtained by substituting equations (B16) into equa- 
tions (B14) and (B15). If only linear terms in the buckling displacements (subscript B 
terms) are retained, and if equations (B17) are utilized, the following buckling equations 
are obtained: 

Nxt-1 + Nxvt-> ~ ® 

^ByX ^B,y 


Nxvt-. + Nv =0 

^B,x yB,y 


N. 


-M- 


yB 


) (B22) 


- Mvx - Mv + -r=r~ + NxWp, _ - Nx„w . 
^B,xx ^B,xy y^BjXy ^B,yy ^ B,xx A 


-%A'^B,yy + 2NxyWB^^y = 0 


with the boundary conditions 






Nxb = 0 

Ub - 

Nxyg = 0 

^^B = 

Mx„ = 0 
^B 

or 

^B,x 

o 

11 

Wb = 

3 given as in equation (! 


> (B23) 




strains and curvatures defined as 
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II 

^^B,x 

Vb,x 



Zb 

II 

'"B,y + 

R 

^xyB 

"'"B,y 

+ ^B,x + 

II 

■^B,xx 


II 

■'^B,yy 


ti 

■2"'B,xy 


(B24) 


It is assumed that there are no additional thermal moments and thermal forces during 
the buckling process so that 


Ntxr = = %Sb = ^Tr^ = ^Txg = = Mjrg = 0 


’B 


■B 


’B 


B 
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NUMERICAL SOLUTIONS 

The solutions to the prebuckling and buckling equations were found by employing 
the finite-difference method. The equations were formed into a system of second-order 
difference equations. As shown in reference 25, this type of system can be easily solved 
by matrix algebra with a modified Gaussian elimination technique. 

In many thermal buckling problems of technical interest, a knowledge of the inter- 
action between thermal and mechanical loads is desired. Hence, for the present problem 
there are two possible choices for a buckling parameter: (1) the magnitude of the applied 
mechanical load (for example, the compressive load or external pressure) for a fixed 
temperature distribution or (2) the magnitude of a temperature change which determines 
thermal forces and moments for fixed mechanical loading. In the present study a solu- 
tion was developed for each parameter. The solution using mechanical loading as the 
buckling parameter is similar in form to that of reference 17 so that the details are 
omitted herein. The solution with temperature as the buckling parameter is described in 
detail in the following sections. 

To employ temperature as a buckling parameter, a scalar temperature change T^ 
is defined so that Tgj^ = T^^/Sgj^, Tg = T^^jSg, and Tp = T^^^; thus, in the shell 

T = TA^sk(l + (Cla) 

whereas in the stringers 

T = TA/3s^l + rse^'^j (Clb) 

and in the rings 

T =Ta^p^ 1 + (Clc) 

Note that in the integrals associated with equations (Bll) and (B12), T^ can be con- 
sidered as a constant. 

By virtue of equations (Cl), the magnitude of the relative temperatures in the stiff- 
ened cylinder is linearly related to a scalar T^. The shape of the temperature distri- 
bution is specified by the constants /3, y, and g. Thus, the problem to be solved is one 
in which the magnitude of the base temperature T^ niay vary while the shape of the 
temperature distribution remains constant. 
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where 


c-od) = — i — dz 

1 - MxMy ^ 

Nt,(1)=^Ct,(i*v8=‘/=) 
with the integration to be performed over the shell wall and 


6iH = 1 

(i = H) 

o 

II 

a 

(i^H) 


The subscript H denotes the difference station at which a ring is located. The barred 
• thermal terms are found by considering a unit temperature change and letting 

"^sk ” ^sk 
Ts = ^s 


Tj. - 


in equations (Bll) and (B12). 

In considering the boundary conditions associated with equation (C4), it was assumed 
that the cylinder was supported at each end (w = 0) and that the ends were free to expand 
longitudinally and a compressive load is prescribed at the boundary. Under these 
assumptions 


Nx^ = -Nx (C5) 

where Nx is the magnitude of the applied compressive load. If the cylinder is supported 
so that w^ = 0 at the ends, the boundary conditions on w (eqs. (B18c) and (B18d)) can 
be stated in finite -difference form at stations i = 0 and i = k (at the ends of the cylin- 
der where x = 0 and x = a, respectively) as follows: For simply supported shells 


^1^0 = % + 


(C6a) 
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and for clamped shells (by using the definition of | and eq. (C3)) 


where 


and 



^6 



fll\ 

CllJ 


^7(i)=§^(NTs + NTxj-(MTx 



(C6b) 


(i = 0,k) 


The governing equations (C4) and boundary conditions (C6) can be written in matrix 
form as 


where 


L K 
I Lj I 

I L2 I 


(Za)o 

(2a)i 

(Za)2 


T^Rl + S 
TaR2 + S 


(C7) 


I Lk_i 
K 



"^A^k-l + S 



a2^2 
2 + 1 

/3l 

-a2 



-2 


(i = 0,1,2,. . .k) 


(i = 1,2,3,. . .k-1) 
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Ri = 



S = 



I = 


1 

0 


0 

1 


(i = 1,2,3,. . .k-1) 


For simple support boundary conditions 


L = I 


K = 


0 

0 


0 

0 



s = 


h 


0 


Rq = 


%( 0 ) 

h 


0 


Rk - 


/37(k) 


0 


and for clamped boundary conditions 



2 

1 
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Because eqioation (C7) contains a band matrix, to conserve computer storage it is 
convenient to solve for the unknowns by Gaussian elimination. By performing a 

series of elementary row operations (see ref. 17), equation (C7) can be written as 



where 


PO = L-% 

Pi = (Li - Pi-l)"^ 


(i = 1,2,. 



(C9) 


and 


Qo = SL“^ =Qk 

Qi = Pi(s - Qi.i) 

(Q]^o “ ^L ^ 

(Qr)i = Pi|Pi ”^i-l] 

(^T)k “ %L ^ 
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The unknowns are determined by finding Pj, Q[, and from equa- 
tions (C9), (CIO), and (Cl 1) and solving successively for (z^jj_2> • • - from 

equations (C8). 


It should be noted from equations (C8) that the prebuckling radial displacements w^ 
are linear with so that 





+ Ta 

1 

("Am), 




(C12) 


where are the variables associated with deformation resulting from mechanical 

loading whereas are the variables associated with deformation resulting from 

thermal loading. Deformations associated with mechanical loading are found by setting 
equal to zero in equation (C8) and solvii^ for Deformations associated with 

thermal loadings are obtained by deleting the Qj column from the right-hand side of 
equation (C8) and solving for (Z^jj corresponding to a unit base temperature change 


Buckling Solution 

If torsional loading is neglected (Njjy = 0^, the partial differential equations gov- 
erning buckling of a stiffened cylinder (eqs. (B22) and (B23)) can be uncoupled into a set 
of ordinary differential equations by assumii^ the following relations: 


Ug = U(x) cos ^ 


Vg = V(x) sin ^ 

) 

Wg = W(x) cos ^ 
Mxg = M(x) cos , 


(C13) 


where n, the number of circumferential buckling waves, is an integer. The substitution 
of equations (C13) into equations (B22) together with the definition of yields four 

ordinary differential equations in u, v, w, and M. To eliminate derivatives of higher 
order than 2, equation (B20) was employed together with the expression 
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X ^11 ^12 1 ^12 ^11/ \ ^ 

"^B,xxx "'Dll ^ ■*■ R ^ ■" " D 


-i^w 
Dll ®»^yy 

(C14) 


obtained from the derivative of Mxg. The resulting buckling equations can be written as 
AllU" + A12U + A13V' + ^Aj 4 - T;^Ai5jw" + ^Aie - TAAi7jw’ 

+ (^18 - TaAi9)w + AiioM’ = 0 


(Cl 5a) 


A21U' + A22V” + A23V + A24W" + ^A25 - TAAggjw' + (a 27 - TaA2s)w = 0 (C15b) 

31 - TaAs2)u' + A33V" + ^A34 - TaA3s^V + ^A 3g - TAA37jw” + (A38 - TAA39 
■ '^1^310)^' + (a3u - TaA312)w - M" = 0 (C15c) 

A41U + A42V + A43W" + ^A^4 - TaA45^W' + A40W + M = 0 (C15d) 


where 


All = ^11 


K 


11 


Dll 

2 


Ai2 = 


Ai3 - ^12 + C 00 - 


Kl2Kii\ n 

Dll R 


Ai4 = ^ll^M 
^15 = -^ll^A^ 

Ai6 = Aiiw;;^+^|Ci2- 


KiiKi2\ 


’ll 


+ [Ki2 + 2K00 - 


^llDl2 yn\^ 

Dll yW 
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Ai7 = 
^18 = 
-^19 = 
^110 
A21 = 
A22 = 
A23 = 
A24 = 
A25 = 
A26 = 
A27 = 
A28 = 
A3I = 
A32 = 
A33 = 





lEil 

Dll 


-(ci2 + Cee)! 



(k,2 + 

-A2l''kT 



'22 


n 


R 



Cfifi - wV 
66 A-r 


■'12 

R 


2 

- (ki 2 t 2K66)(|) -Ciiwk^ 


■ll^A 


T 


-2K 


n 

66 R 
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A34 = 

^35 = 
^36 = 

^37 = 
^38 = 

^39 = 
^310 

^311 

^312 


^22“ . ,, /n 

iR 


R 


2 -Ms) 


C 1 0 — wV 
12 R A-j- 


KiiwX 


K 


12 


R 


M 


+ N, 


GgJg V 

4D06 + ^ + 2^ 6(x - jl) Gj-Jr + D 

j=l 


12 


■*^11^At 

^12 ^Kj2 + 2K66)(|' 


R 

N 


w 


' ...M 




- 2 6'(x - jl) GrJr(| 

j=l 




Ki2 + 2K66)^^ 


R 


w'a + Cii/wV w'. + w' w'l \ 

At Am Axj 


^h'^At^^At 


C12K11 


'll 


2 Kj2 - K66, 





^ll/W \^l R W r2 CuVR/ 


( k ,2 - K 66 ) - - ^^2 


^(sKt - [«Tr + NTy - 



H' 

R/ 
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A41 = -Kii 

A42 = -Ki2 ^ 

A43 =®11 
A44 = 

A45 = KilW^^ 


A46 




To represent the Dirac delta and doublet functions appearing in the coefficients of 
equation (C15), the following finite-difference approximations were employed: 

6(x - jl) f = 6iH ^ 


5’(x - jl) f = 5(i_i)H^ - 6(i+i)H^ 

In this form, equations (C15) are the analogs of equations (36) of reference 17. As 
outlined in that reference, equations of this form can be cast into matrix form, and buck- 
ling loads can be found by examination of a simple determinant resulting from a Gaussian 
elimination procedure similar to that previously described herein. A solution of this type 
was obtained in the present study by using the determinant plotting method described in 
reference 17. 


Modal Iteration Scheme 

If the central difference formula (eqs. (C3)) are employed for the first and second 
derivatives of U, V, W, and M in equations (CIS), the equilibrium equations can be 
written at the interior stations of the shell as 


AiZi_i + BiZi + CjZj^j - 


"^A^^^i^i-l + ®i^i ^i^i+l] 


(i = 1,2,. . .k-1) (C16) 
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where 



All 

A . . A . 

-■g Ai 3 Ai4--Aig 

-f Auo 


Ai - 

-f A21 

A22 A 24 - j A 25 

0 


-f A3I 

A33 A30 - j A38 

-1 



■f ^41 

0 A43 - y A44 

0 


- 2 Aii + a 2 Aj 2 

0 “^-^14 

+ A^Ais 

0 


0 

-2A22 + a2A23 -2A24 

+ a2A27 

0 


0 

-2A33 + A^A 34 -2A38 + a 2 a 3 ii 

2 


0 

A^A42 -2A43 + a2a48 

a 2 


All 

f Ai3 

Ai4 + ^ Ajg 

A A 

^ Alio 

f A21 

A22 

A24 +f A 25 

0 

f A31 

A33 

A36 A38 

- 

■ 1 

|a41 

0 

A43 J A 44 

0 


0 

0 

A15 -f Aj7 

0 



0 

0 

-f A 26 

0 


D. = 






1 

-f A32 

0 

A37 A39 

0 



0 

0 

-f A45 

0 
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- 2 Ai 5 + A^Ajg 

a2a28 

-2A37 + a2a312 
0 

^15 +f ^17 
f A26 

A37 + j A39 

f A45 
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with 


^ /n \2 


^23 = "^ 22 ! 


A27 = -C66 I " ^22 ^ - ^22(11 


R 


A 34 = C 22 ^-^K 22 |) 


R 


n ,„t. 




G«J, 


Gt>Jt 


v 2 K 


A 36 = - ( 4 D 66 + ^ + «iH ^ 


A 38 = 


~R * (‘^12 ■" ^"^6^(1) '"Am - CiiWam*Am * °rJr(|f [6(1-1)H ' «(1+1)H 


A 311 °| ‘^c\^“ - 2 (Ki 2 


K, 


66 


n \2 ^12 



i 77 „ ^ 12 Vn' 


R 




R/ 


^ 2^/^\2 ^ C22 Ci 2 /n \2 


R R 





e) ■ 


C 12^11 


11 


^r] RJ At ■ Rr22 - Cj^^r] "^At 


N»Pj. + Nrpy 


^(nts-n 4 (|- 


in which 


C 22 = 


K 22 = 


D 22 = 




dz + 6iH 


ErA 


r^r 


1 - 


1 - 


Sy E-pA-nZii 

z dz + 6iH 

A 

Ey - /Eyl- E„A_z 5 ^ 

y z2dzH-6iH-p+-£^J 
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In order to develop a modal iteration solution, the following forward and backward 
difference formulas are introduced at ends of the cylinder (finite-difference stations 0 
and k) 




(C17) 


The cylinder was assumed to be supported at the ends so that Wg = 0, The eight possible 
boundary conditions for buckling in equations (B23) can be written in matrix form as 


BqZo + CqZi = T^EqZi 

■^k^k-1 ®k^k = "^A^k^k-l I 


(C18) 


with 


Bq = 0! 


-All 


0 

Dll 


l-Q-ll 

0 

0 

0 

-C66|a 

•■^66 

0 

0 

+ 

0 

1 -■^2 

0 

0 

0 

0 

0 

0 


0 

0 

1 

0 

0 

0 

-1 

0 


0 

0 

0 

1 - -544 


(C19a) 


Cq = a 


A 


11 


'66 




2K66| 


(C19b) 
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where 


0 0 -AiiwX^(O) 0 


0 


Eq =0! 


0 


0 

0 


0 

0 


0 

0 


0 


0 


0 


0 


and 



0 0 0 



«22 ® 

0 0 


0 

0 


0 


0 0 Q!44 


Ak = Co 
Bk = Bo > 
Ek = -Eo^ 


where in equations (C20) w^ terms are to be evaluated at x = a rather than at 
The elements of the selection matrix a are defined as follows: 


Element 

Value 

Boundary condition 

“11 

0 

Ub = 0 


1 

Nxg = 0 

^2 

0 

o 

II 

PQ 

> 


1 

Nxyg - 0 

“44 

0 

o 

II 

1 

1 

o 

II 

> 


(C19c) 


(C20) 


X = 0. 
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The equilibrium equations (C16) and boundary conditions (C18) can now be written 
in matrix form as 


AX = T^BX + T^CX (C21) 

where is a scalar, X is a vector, and A, B, and C are matrices. In order to 
perform modal iteration, let 



so that equation (C2 1) can be written as 


(C22) 


XAX = BX + - CX (C23) 

X 

To start the iteration process assume that X = X° where X° is a trial eigenvector. 
Substitute X° into the right-hand side of the approximate equation 

XAX «BX 


and solve 


AXX = BX° (C24) 

by treating XX as the unknown variable. In finding the solution to equation (C24), 
Gaussian elimination (Potter’s method) is employed to avoid having to find the inverse 
of A. The resulting approximate eigenvalue Xj and normalized eigenvector X^ can 
be taken as 

Xi = ixxl 

xl = — (C25) 


For subsequent iterations the equation 

AXX = BXi + — CXi (i = 1,2,. . .) (C26) 

H 

must be solved to find more accurate solutions to equation (C23). The iteration process 
converges to the largest characteristic value of X (lowest T^) provided - CX is 

A 

small compared with BX. Thus the temperature T^ can be found for a specified num- 
ber of circumferential waves n. The lowest value of T^ found is then the buckling 
temperature. 
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COMPUTER PROGRAMS 
BAMSOC I 

The computer program BAMSOC I finds the buckling load for a stiffened cylindrical 
shell subjected to mechanical and thermal loadings by using determinant plotting. The 
programed solution is similar to that developed in reference 17 and differs only in that 
the present solution considers thermal effects in a multilayered orthotropic cylinder. 
Buckling loads or buckling temperatures can be selected as the buckling parameter, and 
critical values are found by determinant plotting over user-specified ranges of circum- 
ferential wave number n and the buckling parameter. For a specified n, the program 
seeks a change in sign in a determinant or "modified residual" by varying the buckling 
parameter. A flow diagram for BAMSOC I (determinant plotting solution) follows: 
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The program is written in FORTRAN IV and will run on the Control Data 6000 
series computer with the Scope 3.0 operating system. The program elements were over- 
layed so that a maximum of 200 finite-difference intervals could be used on a remote 
computer terminal with 70 000 octal storage units. A maximum of 20 layers can be used 
in the shell. Input required by the program is contained in the following namelists: 


$ CON 
MI 

NF 

DELN 
NBARXI 
NBARXF 
DELNX 
PINIT 
PFIN 
DELP 
TAI 
TAF 
DELTA 
KASE = 1 
= 2 
= 3 
= 4 

$ (end of CON) 
50 


initial value of n to be investigated 
final value of n 
increment in n 

A 

initial or fixed value of N^, axial compressive load 
final value of N^ 

A 

increment in Nx 

initial or fixed value of p or q, lateral or hydrostatic pressure 
final value of p or q 
increment in p or q 

initial or fixed value of T^, base temperature (see eq. (Cl)) 
final value of T^^ 
increment in T^^ 

axial compression with fixed temperature and lateral pressure 
lateral pressure with fixed temperature and axial compression 
hydrostatic pressure with fixed temperature and axial compression 
thermal loads corresponding to T^ with fixed mechanical loads 


II 
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$ NANCY 
R 
A 

DELBAR 
BC = 1 
= 2 
ER 
AR 
IR 
ZR 
GRJR 
NRING 
NRSPACE 

ES 

AS 

IS 

zs 

GSJS 

D 


radius to reference surface 
cylinder length 

distance from inner wall surface to reference surface 
simply supported for prebuckling analysis 
clamped for prebuckling analysis 
Young's modulus for ring 
area of ring 

moment of inertia of ring about its centroid 

distance of centroid of ring to reference surface (positive, external rings) 
ring torsional stiffness, GJ 
number of rings 

number of finite-difference intervals between adjacent rings (no more 
than 200 stations along cylinder length) 

Young's modulus for stringer 

area of stringer 

moment of inertia of stringer about its centroid 

distance from centroid of stringer to reference surface (positive, 
external stringers) 

stringer torsional stiffness, GJ 

stringer spacing 


51 



APPENDIX D - Continued 


ECC 

distance from reference surface to surface on which Nx 3^cts 

GAMS 

stringer temperature shape parameter, yg (see eq. (Clb)) 

GAME 

ring temperature shape parameter, (see eq. (Clc)) 

TAUS 

y o^^s(^)dAs 


TAUR 

TAUSPR 

^ Q!j./3j.(z)dAr 

y ois)3s(z)z dAs 

where jS is temperature shape parameter in 
^ equation (Cl) 

TAURPR 

J ari3r(z)z dAj. 


G 

exponential temperature shape factor, g (see eq. (Cl)) 

EX(I) 

Ex of ith layer of shell (first layer lies on inner surface of shell 
more than 20 layers) 

EY(I) 

Ey 


NUX(I) 

Mx 


NUY(I) 

Py 


GXY(I) 

^xy 


H(I) 

thickness of layer 

GAMMA(I) 

layer temperature shape factor ygj^ (see eq. (Cla)) 

BETA(I) 

layer temperature shape factor jSgj^ (see eq. (Cla)) 

ALPHAX(I) 

«x 


ALPHAY(I) 



LAYER 

number of layers (^0) 


$ (end of NANCY) 
52 


APPENDIX D - Continued 


$ CORNIE 


lALF = 1 

Nxb = 0, Vg = 0, Wg = 0, Mxg = 0 

= 2 

Ur ~ ''^B ” ^ 

= 3 

Nxb = 0, Nxy^ = 0, Wg = 0, Mxg = 0 

= 4 

Ug = 0, Nxyg = 0, Wg = 0, Mxg = 0 

= 5 

Ug = 0, Nxyg = 0, Wg = 0, Wg = 0 

= 6 

Nxg = 0, Nxyg = 0, Wg = 0, Wg = 0 

= 7 

Nxg = 0, Vg = 0, Wg = 0, Wg = 0 

= 8 

Ub = 0» Vg = 0, Wg = 0, Wg^ = 0 


MODE = 0 calculation of buckling mode shape is not required 

■* 0 calculation of buckling mode shape is required 

IRING = 1 simplified ring theory described in appendixes A, B, and C 
= 2 more accurate ring theory described in appendix E 

$ (end of CORNIE) 

$ RING (required only if IRING = 2) 

IX Ix, moment of inertia for out-of-plane bending of ring (eq. (Ell)) 

lOXZ Ioxz> cross product of inertia for out-of-plane bending of ring (eq. (Ell)) 

$ (end of RING) 

A sample output listing for BAMSOC I follows. The case shown is the ring- and 
stringer-stiffened cylinder loaded by thermal loading alone. (See text.) The parameter 
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has been selected as the change in temperature of the heated skin. Experience with 
determinant plotting solutions has shown that it is possible to miss sign changes in the 
determinant by taking too large an increment in the load parameter. As a safeguard 
against this possibility, the present program examines small loading increments, when 
unusual trends in the determinant are detected. 

In the sample problem, for N = 51 or 52, the modified residual value was found to 
increase in value after a steady decline. The irregularity was detected and finer load 
increments were employed in the calculations and a change in the sign of the determinant 
was found. Note, however, that no irregularities were found for N = 53, and as a con- 
sequence, sign changes were found at higher temperatures. Further investigation with 
finer load increments in calculations not presented established a buckling temperature of 
about 525° F (292 K) for N = 57. 
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namelist input for BAMSOC I 


SCON NI*50*» NF*53*« OELN= 1 • ♦ NBARXI=0** NBARXF=0** OELNX*!** 

PINIT«0#* PFINaO** DELP=1*, TAI*450*« TAF=600*» DELTA=ao*» KA5E=4S 

SNANCY R=65.f A«300.. DELBaR“.03* 0C=1« ER»16.4E6* AR= • 246 ♦ ‘ I R* • I 377 1 
ZR*-1#992* GRJR=3477«* NR1nG*19* NRSPACE=10* ES=16*4E6» AS**C495« IS=*004082* 

ZS*-*342« GSJSa93*« D=l*35« ECC=0*» GAMS»0*« GAMR=0* * TAUS=0*i TAUR=0#< 
TAUSPRaO** TAURPRaO** G=0#» EX ( 1 ) = ? 4 .SEB ♦ EY < 1 ) = I 4 . 5E6 « NUX ( I ) = • 3 • NUYtl) = .3* 
GXYd )«5*580E6t H(l)=#06* GAMMA(l)aO#« BETA(1)=1«* 

ALPHAXU )a*5.E-6* ALPHAYd )=S.E~6* LAYER=IS 
SCORNIE IALF=1 ♦ MODE*©* lRiNG=t * 


BAMSOC I 


CALCULATION OF BUCKLING LOAD FOR THERMALLY STRESSED. STIFFENED* MULTI-LAYERED CYLINDER USING DETERMINANT PLOTTING 
CHANG-CARD R0F364 NASA LANGLEY RESEARCH CENTER 1968 

DATE 03/27/69 


NI= 5.0000E+01 
TAI* 4.5000E+02 
NBARXa 0. 

P= 0. 


NF= 5.3000E-I-01 0ELN= 1 .OOOOE+00 

TAF= 6,0000E+02 DELTA* 2.0000E+01 


INPUT FOR PROGRAM THERMUL 

CHANG-CARD 

A2031 R0F364 JANUARY, 1968 



GEOMETRY PARAMETERS 

LOADINGS 

TEMPERATURE PARAMETERS 



R* 6.5000E+01 


TAUS= 0. 

GAMS= 

0 

A= 3.0000E-»-02 


TAUR= 0. 

GAMR= 

0 

DELBAR* 3.0000E-02 

ECC= 0. 

TAUSPR= 0. 

G= 

0 



TAURPR* 0. 




RING PROPERTIES 

ER= 1.6400E+07 ZR=-l,9920E+00 
AR= 2.4600E-01 GRJR* 3.4770E*03 
IR= 1.3770E-01 


STRINGER PROPERTIES 
ES= 1.6400E<-07 
AS= 4.9500E-02 
IS= 4.0820E-03 


ZS=-3.4200E-01 
GSJS= 9.3000E+01 
D= l,3500E+00 


NUMBER OF RINGS* 19 

NUMBER OF FINITE DIFFERENCE SPACES BETWEEN RINGS* 10 


GEOMETRY AND TEMPERATURE PROPERTIES FOR LAYERS 


LAYER EX EY NUX 

1 1.450E+07 1.450E+07 3.000E-01 

SIMPLE SUPPORT BOUNDARY CONDITIONS 
11 12 


NUY 

3.000E-01 

GXY 

5.580E+06 

H 

6.000E-02 0 

22 


66 


GAMMA 

BETA 

ALPHA X 

ALPHA Y 


l.OOOE+00 

5.000E-06 

5.000E- 


13 23 


B 1 l,59340659E+07 

C 9.56043956E+05 

0 2.86813187E+02 

K 0. 

L 

M 


4.78021978E•^06 

2.86813187E+05 

8.60439560E4-01 

0 . 


1 , 59340659E+07 
9.56043956E+05 
2. 86813187E+02 

0 . 


5.58000000E+06 
3.34800000E+05 
1 .004^0000E+02 

0 . 


-1 .03571429E+02 -1 . 0357 1429E + 02 


-6.21428571E^•00 -6 . 2 1 4285 71 E + 00 

0 . 0 . 
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INPUT FOR CALCULATION OF 
IALF= 1 


N* 5.0000E+01 
TA 

4.50000000E+02 
4.70000000E+02 
4. 90000000 E+02 
5.10000000E+02 
5.30000000E4-02 
5.50000000E+02 
5.40000000E+02 
5.40000000E+02 


FIEO RESIDUAL 


MODIFIED RESIDUAL 

-7.79225885E>05 
-7.75252974E«-05 
-7.6S576851E+05 
-7.6057l046E<-05 
-7.42188514E+05 
6.79009921E<-05 
-7.14898570E<-05 
CROSSING FOUND WITHIN 1 


CHANG-CARD RDF364 


PERCENT 


N= 5.1000E+01 


TA 


MODIFIED RESIDUAL 


4.50000000E+02 
4.70000000E+02 
4.90000000E+02 
5.10000000E^02 
5.30000000E+02 
5.50000000E+02 
OPEN UP INTERVAL IN 
5.30Q00000E+02 
5.32000000E+02 
5.34C00000E+02 
5.36000000E+02 
5.38000000E+02 
5.40COOOOOE+02 
5.39000000E+02 
5.39000000E+02 


-7.88803726E<-05 
-7.84474819E+05 
-7.78172671E+05 
-7.677868 96E4-0 5 
-7.43541985E*05 
-9.24281493E+05 
SEARCH CF CROSSING 

-7.43541985E+05 
-7,38292392E^05 
-7.31168334EI-05 
-7.19848153EI-05 
-6.91640017E+05 
9.03024088E+05 
-6.24795117E+05 
CROSSING FOUND WITHIN 


PERCENT 


N= 5.20C0E4-01 


TA 


MODIFIED RESIDUAL 


OP 


4.50000000E+02 
4•70000000E^-02 
4.90000000E+02 
5.10000000E+02 
5.30000000E+02 
5.50000000E+02 
EN UP INTERVAL IN SEARCH 
5.30COOOOOE+02 
5,32000000E+02 
5.34000000E+02 
5.36000000E+02 
5.35000000E+02 
5.35000000E+02 Cl 


-7.98256606E^-05 
-7. 9359432 4E <-05 
>'7.866824O0E<-O5 
-7.74838733Ef05 
-7.41093951E<-05 
-7.616 22322E+05 
}F CROSSING 
-7.4l093951E<-05 
-7.29620508E<-05 
-6.99861698E<-05 
8.61l74311E<-05 
-6.16100274E<-05 
ISSING FOUND WITHIN 1. 


PERCENT 


N= 5.3000E+01 


TA 


MODIFIED RESIDUAL 


4.50000000E<-02 

4.70000000E+02 

4.9C000000E+02 

5.10C00000E+02 

5.30000000E+02 

5.50000000E+02 

5.70000000E<-02 

5.60000000E+02 

5.60000000E<-02 


-8.07594308E<-05 

-8.02630842E+05 

-7.95147042E<-05 

-7.818l7375E<-05 

-7.24361857E»05 

-7.05664617Ef05 

6.79333836E+05 

8-50808776Ei'05 

CROSSING FOUND WITHIN 1. PERCENT 


MINIMUM TA= 5.35000000E<-02 


N= 5.20000000E<-01 


FEBRUARY, 1968 
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APPENDIX D - Continued 
BAMSCX: n 

The computer program BAMSOC n finds the buckling load for a stiffened cylindrical 
shell subjected to mechanical and thermal loadings by using modal iteration. The solu- 
tion corresponds to equations (C16), (C18), and (C27) of the present text. A flow diagram 
for BAMSOC n (modal iteration solution) follows: 



Buckling temperatures are found for specified ranges of the circumferential wave num- 
ber n. If the solution is well behaved, the program will converge to the candidate buck- 
ling temperature after several iterations. The program is overlayed to permit a maxi- 
mum of 200 finite-difference intervals and requires 150 000 octal storage units. 
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APPENDIX D — Continued 


The input required by the modal iteration solution is similar to that for BAMSOC I 
except that loading ranges and increments are not specified and the solution is valid only 
for simply supported boundary conditions. The following modifications are required for 
the three namelists presented for BAMSOC I: 

(1) in namelist CON, retain only NI, NF, DELN, NBARXI, and PINIT 

(2) retain namelist NANCY 

(3) in namelist CORNIE, permissible values of lALF are from 1 to 4. Delete 

MODE and IRING. 

The convergence criteria (see eq. (C27)) is specified internally for the ith iteration as 


K^'a) 

i - (Ta)i.1 


(Ta)i 



0.001 


A sample output listing for BAMSOC n follows. The case is the same case as for 
BAMSOC I except that only 100 finite-difference intervals were employed to condense the 
output. 
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APPENDIX D - Continued 


NA'ltLlST INPUT FOR BAMSOC H 


SCON NI=57.* NFsST,* 0ELN=1.* NdARXl=0.* PlNlT=0.S 

SNA.NCY P=6b.t As300,* OtLdAP=,'Jj* dC*lt ER=16,AE6* ARs.d46f IR=.l377f 
ZR=-l,992t GRJRS347/.. NRING=19* NRSPACE= 5* ES=lto*4£6* AS=.0495, IS=.004082» 
ZS=-.342f GSJS=^3.* 0=1.35* LCC=0.* GAmS=0.* GAMR=0.* TAUS=0.* TAuR=0.* 

TAUSPR = n,, TAUR3R = 0.* G=0.* GXY ( 1 ) =b . 5ttOE6* H(l) = .06* GAHJIA ( 1 ) =0 • , BETAU)=l.* 
EX<1 )s 1A,SE6, EY(l)=l4.bE6* NuA(l>=.30* NUY<l)=.3u* ALPHAX ( U -5.E-6, 
ALPHAY(l)=5.t-6. LAVtK=lS 

5CORNIE IAlFsIs 


BAHSOC II 

CALCULATICN OF BUCKLING LOAD FOR THERMALLY STRESSED, STIFFENED, MULTI-LAYER CYLINDER BY MODAL ITERATION 
CHANG-CARO ROF364 NASA LANGLEY RESEARCH CENTER 1968 

DATE 03/27/69 

NI» 5.7000E+01 NF« 5.7C00E+01 DELN» l.OOOOE+00 


INITIAL Z VECTOR (UNSCALEOl 


UB 


V6 


MB 


MB 


1 0 . 0 . 

21 0 . 0 . 

41 0. 0. 

61 0. 0. 

81 0. 0. 

101 0 . 0 . 


l.OOOOOOOOE+OO 0. 
1.00000000E«'00 0. 
1.00000000E4*00 0. 
1.00000000E<-00 0. 
l.OOOOOOOOE+00 0. 
l.OOOOOOOOE^OO 0. 


ITERATION NUMBER 1 


TEMPERATURE TA* 1 . 63659058E-C 1 CORRESPONDING SCALED Z VECTOR 



UB 

VB 

WB 

MB 

1 

-7.76570953E-07 

0. 

0. 

0. 

21 

-1.91130069E-10 

1. 83146093 E-08 

7.35752664E-09 

-1.83573040E-01 

41 

-7.01489655E-14 

1.83760317E-08 

7.39760844E-09 

-1 .83772776E-01 

61 

7.01489688E-14 

1.83760317E-08 

7. 39760844E-09 

-1 .83772776E-01 

81 

1.91130C69E-10 

1.83146093E-06 

7.35752664E-09 

-1.83573040E-01 

101 

7.7657C953E-07 

0. 

0. 

0. 


ITERATION NUMBER 2 


TEMPERATURE TA» 1 . 1 7890697E<'0 3 CORRESPONDING SCALED Z VECTD^ 



UB 

VB 

W6 

MB 

1 

7.27927938E-07 

0. 

0. 

0. 

21 

-1.57725826E-09 

5.04390724E-08 

3.52955211E-08 

-1. 52729 564E-01 

41 

-1.2039 8426E-12 

5.07973257E-08 

3.55277678E-08 

-1 .53943644 E-01 

61 

1.20398423E-12 

5.07973257E-08 

3. 55277678E-08 

-1.53943644E-01 

81 

1.57725826E-09 

5.04390724E-08 

3.52955211E-08 

-1.52729564E-01 

101 

-7.27927938E-07 

0. 

0. 

0. 


ITERATION NUMBER 25 


TEMPERATURE TA- 5. 37719505E*02 CORRESPONDING SCALED Z VECTOR 



UB 

VB 

WB 

MB 

1 

4.58064616E-07 

0. 

0. 

0. 

21 

-2.87625169E-07 

-9.81722985E-09 

-6.56022956E-09 

3.23334452E-02 

41 

-4. 84982682 E-08 

-4.40150651E-09 

-2.92732002E-09 

1.45714341E-02 

61 

4.04982682E-O8 

-4.40150651E-09 

-2.92732002E-09 

1.45714341E-02 

81 

2.87625169E-07 

-9.81722985E-09 

-6. 56022956E-09 

3.23334452E-02 

101 

-4.58064816E-07 

0. 

0. 

0. 



APPENDIX D - Continued 


ITERATION NUMBER 26 


TEMPERATURE TA= 5.37123003E+02 CORRESPONDING SCALED Z VECTOR 



UB 

VB 

HB 

MB 

1 

4.46578687E-07 

0. 

0. 

0 • 

21 

-2.93047826E-07 

-9.38196491E-09 

-6.26962170E-09 

3. 08988773 E-02 

41 

-5.29982857E-08 

-4.69946566E*09 

-3.12613561E-09 

1.55542013E-02 

61 

5,29982857E“08 

-4.69946566E-09 

>-3.12613561E-09 

1.55542013E-02 

81 

2.93047826E-07 

-9.38196491E~09 

-6.26962170E-09 

3.08988773E-02 

101 

-4.46578687E-07 

0. 

0. 

0. 


ITERATION NUMBER 27 


TEMPERATURE TA= 5.3657592AE^^02 CORRESPONDING SCALED Z VECTOR 



UB 

VB 

WB 

HB 

1 

4.35668C83E-07 

0. 

0. 

0. 

21 

-2.97920426E-07 

-8.94910655E-09 

-5.98057394E-09 

2.947250 94E-02 

41 

-5.75110986E-08 

-4.98896109E-09 

-3.31932851E-09 

1.65089071E-02 

61 

5.75110986E-08 

-4.98896109E-09 

'3.31932851E-09 

1.65089071E-02 

81 

2.97920426E-07 

♦8.94910655E-09 

-5.98057394E-09 

2.94725094E-02 

101 

-4.35668C83E-07 

0. 

0. 

0. 


ITERATION NUMBER 28 


CONVERGED TEMPERATURE TA= 5. 3607A232E+02 CORRESPONDING SCALED Z VECTOR 



UB 

VB 

HB 

MB 

1 

4.25301149E-07 

0. 

0 . 

0. 

2 

-7.59850632E-07 

-4. 59 62 5262 E-07 

2.88703579E-05 

1.97151230E-01 

3 

-5.21414819E-07 

-7.71134756E-07 

4.27008188E-05 

2.03472475E-01 

4 

7.46228255E-11 

-7.28756070E-07 

3.92369407E-05 

1.45920054E-01 

5 

2.70353608E-07 

-4.414841706-07 

2.31230666E-05 

8.13483855E-02 

6 

2.99574280E-07 

-9.37078924E-09 

-5.52645553E-09 

-6.63765673E-04 

7 

2.81086829b-07 

4.28726627E-07 

-2.32819862E-05 

-1.08527967E-01 

8 

1.11722331E-07 

6.97014670E-07 

-3.79640995E-05 

-1.76416032E-01 

9 

-1.45090817E-07 

7.09547922E-07 

-3.85226372E-05 

-1.84431761E-01 

10 

-3.27532512E-07 

4.55146445E-07 

-2.42795296E-05 

-1.24870198E-01 

11 

-3.62973005E-07 

2.43286648E-09 

1.72226769E-09 

-9.63841827E-03 

12 

-3.44911857E-07 

-4.67981720E-07 

2.50584371E-05 

1.20262254E-01 

13 

-1.44970992E-07 

-7.46556218E-07 

4.04388090E-05 

1.92473613E-01 

14 

1.49747724E-07 

-7.43563747E-07 

4.02956732E-05 

1.90551988E-01 

15 

3.45113853E-07 

-4.61421596E-07 

2.47843487E-05 

1.15363571E-01 

16 

3.53028614E-07 

4.79054708E-09 

3.21738900E-09 

-1.57017051E-02 

17 

3.03512382E-07 

4.42124012E-07 

-2.34938675E-05 

-1.23089686E-01 

18 

1. 11311952 E-07 

6.84191543E-07 

-3.70558546E-05 

-1.78314502E-O1 

19 

-1.52077748E-07 

6.67749535E-07 

-3.62463455E-05 

-1.68109183E-01 

20 

-3.16278412E-07 

4.04730815E-07 

-2.188995546-05 

-9.40429520E-02 

21 

-3.02289280E-07 

-8.52149244E-09 

-5.69499266E-09 

2.80636251E-02 

22 

-2.40603980E-07 

-3.70626319E-07 

1.958333956-05 

1.08943509E-01 

23 

-7.79210795E-08 

-5. 60691973 E-07 

3.D3633453E-05 

1.47419054E-01 

24 

1.33393533E->07 

-5.38503291E-07 

2.92702575E-05 

1.33661345E-01 

25 

2.58710335E-07 

-3.20106295E-07 

1.74149936E-05 

6.96391649E-02 

26 

2.34182913E-07 

9.58375632E-09 

6.39701603E-09 

-3.159688036-02 

27 

1.73529469E-07 

2.816430896-07 

-1.48179466E-05 

-8.67100390E-02 

26 

4.90093406E-08 

4.17721507E-07 

-2.26210193E-05 

-1.10635773E-01 

29 

-1.05031425E-07 

3.953797576-07 

-2.15199960E-05 

-9. 67952493 E-02 

30 

-1.92160659E-07 

2.30933170E-07 

-1.26321987E-05 

-4.70731332E-02 

31 

-1.65745567E-07 

-8.73802977E-09 

-5.82660431E-09 

2. 884280696-02 

32 

-1.14405427E-07 

-1.96104461E-07 

1.02615469E-05 

6.28887408E-02 

33 

-2.73157633E-08 

-2.849283886-07 

1.54304282E-05 

7. 599497 09E-02 

34 

7.54101910E-08 

-2. 658063066-07 

1.44877480E-05 

6.41582849E-02 

35 

1.30679669E-07 

-1.52498812E-07 

8. 388828296-06 

2.89144603E-02 

36 

1.07290082E-07 

7.01425061E-09 

4.67248645E-09 

-2.31792101E-02 

37 

6.81917701E-08 

1.24588002E-07 

-6. A8369788E-06 

-4.18076618E-02 

38 

1.24C77770E-08 

1.76889287E-07 

-9.57824365E-06 

-4.76148396E-02 

39 

-4.98547675E-08 

1.62C46387E-07 

-8.84625794E-06 

-3.84342699E-02 

40 

-8.10843369E-08 

9,071498506-08 

-5.02850473E-06 

-1.53969322E-02 
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APPENDIX D - Concluded 


41 

-6.20122302E-08 

-5.26917618E-09 

-3.50634900E-09 

1.74328862E-U2 

42 

-3.40867935E-08 

-6,97854476E-08 

3.59568563E-06 

2.52369502E-02 

43 

-2.04212480E-09 

-9,505316746-08 

5.14159204E-06 

2.61273885E-02 

44 

3.04004734E-08 

-8.38477792E-08 

4.58879044E-06 

1,920219076-02 

45 

4,40193 554E-08 

-4-4194 1266E-08 

2.49591889E-06 

5,26846401 E-03 

46 

2.79055634E-08 

4-04701034E-09 

2.69043138E-09 

-1. 34039223 E-02 

47 

7.88264544E-09 

2.81461217E-08 

-1.39425150E-06 

-1.29366908E-02 

48 

-6,66980972E-09 

3.22525792E-08 

-1.73154237E-06 

-9.83661853E-03 

49 

-1.66255869E-08 

2.30691997E-08 

-1.27836897E-06 

-4.16443057E-03 

50 

-1.61454443E-08 

7.16C82486E-09 

-4.92203515E-07 

3.44703134E-03 

51 

-5.11607606E-20 

-3.61446716E-09 

-2.40179628E-09 

1.19772649E-02 

52 

1.61454443E-08 

7-16082486E-09 

-4.92203515E-07 

3.44703134E-03 

53 

1.66255869E-08 

2.3069 1997E-08 

-1.27836897E-06 

-4.16443057E-03 

54 

6.66580972E-09 

3.22 525792 E-08 

-1.73154237E-06 

-9.83661853E-03 

55 

-7,80264544E-O9 

2.81461217E-08 

-1.39425150E-06 

-1.29366908E-02 

56 

-2,79055634E-08 

4.04701034E-09 

2.69043138E-09 

-1.34039223E-02 

57 

-4.40193554E-08 

-4.41941266E-08 

2.49591889E-06 

5,268464016-03 

58 

-3.04004734E-08 

-8.38477792E-08 

4,58879044E-06 

1.92021907E-02 

59 

2.04212480E-09 

-9;50531674E-08 

5.14159204E-06 

2.61273885E-02 

60 

3.40867935E-08 

-6.97854476E-08 

3.S9568563E-06 

2.52369502E-02 

61 

6.20I22302E-08 

-5.26917618E-09 

-3.50634900E-09 

1.74328862 E-02 

62 

8.10843369E-08 

9-07149850E-08 

-5.02850473E-06 

-1.53969322E-02 

63 

4.98547675E-08 

1.62046387E-07 

-8.84625794E-06 

-3.84342699E-02 

64 

-l*24077770E-08 

1.766892876-07 

-9.57824365E-06 

-4.76148396E-02 

65 

-6.81917701E-08 

1.24588002E-07 

-6.48369788E-06 

-4.18076618E-02 

66 

“l,07290C82E-r07 

7.01425061E-09 

4.67248845E-09 

-2.31792101E-02 

67 

-l,30679669E-07 

-1.52498812E-07 

8.38882829E-06 

2.89144603E-02 

68 

-7.54101910E-08 

-2.658063C6E-07 

1.44877480E-05 

6.41582849E-02 

69 

2.73157633E-08 

-2.84928388E-07 

1.54304282E-05 

7.59949709E-02 

70 

1.14405427E-07 

-1.961044616-07 

1.02615469E-05 

6.28887408E-02 

71 

1.65745567E-07 

-8.738029776-09 

-5.82660431E-09 

2.88428069E-02 

72 

1.92160659E-07 

2.309331706-07 

-1.26321987E-05 

-4.70731332E-02 

73 

l,05031425E-07 

3.953797576-07 

-2.15199960E-05 

-9.67952493E-02 

74 

-4,9OO934O6E-O0 

4,177215076-07 

-2.26210193E-05 

-1.10635773E-01 

75 

-1.73529469E-07 

2.81 £430896-07 

-1. 48179466 E-05 

-8.67100390E-02 

76 

-2.34182913E-07 

9.583756326-09 

6.39701603E-09 

-3. 15968803 E-02 

77 

~2.5071O335E-O7 

-3.201062956-07 

1.74149936E-05 

6.96391649E-02 

78 

-1.33393533E-07 

-5.385032916-07 

2.92702575E-05 

1.33661345E-01 

79 

7.79210795E-08 

-5.606919736-07 

3.03633453E-05 

1.47419054E-01 

80 

2.40603980E-07 

-3.706263196-07 

1.95833395E-05 

1 .08943509E-01 

81 

3.02289280E-07 

-8.521492446-09 

-5.69499266E-09 

2.80636251E-02 

82 

3.16278412E-07 

4.047308156-07 

-2.18899554E-05 

-9.40429520E-02 

83 

l,52C77748E-07 

6.677495356-07 

-3.62463455E-05 

-1.68109183E-01 

84 

-1.11311952E-07 

6.841915436-07 

-3.70558546E-05 

-1.78314582E-01 

85 

-3.03512382E-07 

4.421240126-07 

-2.34938675E-05 

-1.23089686E-01 

86 

-3.53028614E-07 

4,790547086-09 

3.217309OOE-O9 

-1.57017051E-02 

87 

-3.45113853E-07 

-4.614215966-07 

2.47843487E-05 

1.15363571E-01 

88 

-1.49747724E-07 

-7.435637476-07 

4.02956732E-05 

1.90551968E-01 

89 

1.44970992E-07 

-7.465562186-07 

4.04388090E-05 

1.92473613E-01 

90 

3.44911857E-07 

-4.679817206-07 

2.50584371E-05 

1.20262254E-01 


3,62973005E-07 

2.432866486-09 

1.72226769E-09 

-9.63041827E-O3 

92 

3.27532612E-07 

4.551464456-07 

-2.42795296E-05 

-1. 24870198 E-01 

93 

1.450908l7E'-07 

7.095479226-07 

-3.85226372E-05 

-1.84431761E-01 

94 

-1.1I722331E-07 

6.970146706-07 

-3.796409956-05 

-1.76416032E-01 

95 

-2,ai086829E-07 

4.287266276-07 

-2.328198626-05 

-1.08527967E-01 

96 

-2.99574280E-07 

-9.370789246-09 

-5.52645553E-09 

-6.63765673E-04 

97 

-2.70353608E-07 

-4.414841706-07 

2.31230666E-05 

8.13483855E-02 

98 

-7.46228255E-11 

-7.287560706-07 

3.92369407E-05 

1.45920054E-01 

99 

5.21414819E-07 

-7,711347566-07 

4.27008188E-05 

2.0347247SE-01 

100 

7.59850632E-07 

-4,596252626-07 

2.88703579E-05 

1. 971512306-01 

101 

-4,25301l49E-07 

0. 

0. 

0. 


NUMBER OF ITERATIONS-^ 28 


MINIMUM TA- 5.36074232E*02 N=« 5. 70000000 E«-0 1 

RUNNING TIME APPROXIMATELY 72 SECONDS 



APPENDIX E 


DEVELOPMENT OF MORE ACCURATE RING THEORY 


Near the completion of this investigation. Dr. David Bushnell of the Lockheed 
Missiles and Space Company verified the results of some of the present calculations by 
using a newly developed program for shells of revolution (ref. 26). In doii^ so, he 
obtained results which indicated that the therm al-bucklii^ behavior of the rir^-stiffened 
shells investigated herein could be influenced by the type of ring theory employed in the 
analysis. Preliminary results suggested that the out-of-plane ring bending stiffness 
during buckling was especially important. 

The purpose of this appendix is to present an extension of the present analysis to 
include the first-order effects of out-of-plane deformations in the rii^s during buckling. 
The ring theory is developed by a method proposed in reference 27 to obtain a consistent 
theory for sturdy rings attached firmly to a shell. A key assumption in the development 
is that compatibility of the ring and shell deformations is enforced along a normal 
passing through the centroid of the ring. Thus, local deformations of the ring elements 
are ignored and the ring behaves as a rigid beam. More exact models for rings which 
enforce compatibility only at the point of attachment have been proposed. (See, for 
example, ref. 28.) However, these models lead to extremely complex, unwieldy equations. 


Ring Displacements 


The ring center of shear and centroid are assumed to coincide. The ring displace- 
ments u|., v^, and w^ are specified in local ring coordinates x, y, and z, which are 

parallel to the coordinates of figure 1 but have their origin at the ring centroid. The ring 
displacements u^ , v^ , and w^ are expressed (ref. 28) as 


ut = Uc + Z^j 

vt = Vc + z^2 - ^c,y> 


(El) 


wt = Wc 
V, 


- X)3j 


with /3j = -W(, and ^2 = y ^ where the subscript c denotes the ring centroid. 
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APPENDIX E - Continued 


The ring coordinates y 


and 


z are 



related to the shell coordinates y and z by 

■>v 



(E2) 


Z = Z - Zj. J 

The displacements in the shell in the vicinity of the ring are assumed to be 



To develop a ring theory for sturdy stiffeners, compatibility of shell and ring defor- 
mation is enforced along a shell normal passing through the centroid of the ring. Matching 


ring and shell displacements along x = 0 


^that is, 




, and so forth\ yields 


x=0 


and 


Uc = u - ZfW^x *1 






wt = w + XW V 
r ,x 


(E4) 


(E5) 


Strain Energy 

The strain energy in the N interior rings on the shell is approximated as 
N 2 R 2 

" I Z ^0 "" •CD' ^ Jr^l,y - Er^rTeJ. 6(x - jZ)dx dy 


(E6) 
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APPENDIX E - Continued 


with 




(E7) 


where the introduction of the Dirac delta function permits the replacement of x by x 
in the derivatives appearing in equation (E5). 


Buckling Strain Energy 

By following the procedures outlined in appendix B, the total potential energy of the 
shell, its stiffening elements, and the external loading can be formed, and the principle of 
minimum potential energy (611 = 0) can be used to obtain nonlinear equilibrium equations. 
The equations governing prebuckling and buckling can be derived by separating displace- 
ments into axisymmetric and small asymmetric parts (for example, u = ua + ub). For 
the more accurate ring theory, the prebuckling equations obtained are the same as those 
developed as equation (1) or equation (B21) in appendix B. Thus, the new terms of the 
more accurate ring theory appear only in the buckling equations. 

In order to make a comparison of the present theory with existing ring theories, the 
strain energy during buckling can be written as 

N 

i=i 

where 


- 4 -^ 


'B,xyy 
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APPENDIX E — Continued 


are employed. The integration over the ring area can be accomplished with the following 
definitions: 


I 

y 

I 


E]|* dA]^ — Ej^Aj* 


EfZ dAr = ZyErAr 


E]|»z^dAjr* — ®r^oz 



dAr “ ^'r^oxz 


^ Er^ 

J Erx2dAr = Eyl^j 


xr = O) 

> (Ell) 




where the subscript o is used to denote that the inertias a^id lojjg; com- 

puted about the shell reference surface, z = 0. 

After integration, equation (E8) can be written in matrix form as 


in which 


and 


HrB = 





dx 



Gil 

Gi2 

Gi3 

Gi4 


G22 

G23 

G24 



G33 

G34 

Symmetric 



G44 


(E12) 


(E13) 


(E14) 
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APPENDIX E - Continued 


with 


N 


Gil = X 
j=l 


N 


Gi 2 = 2 / ^ 

j=l 


N 


Gi 3 = ^ 6 (x - J^)’?ErIoxzfe) 
j=l 


N 

Gi4 = ^ 6(x - ]7) 

j=l 

N 

G 22 = X 

j=l 

N 

G 23 " ^ 
j=l 


2 _ / \4 

R ' ^ 


/n \2 - n 2 n 2 

ErAr(— ) + 2 zj.ErAr — - + Ej-Ioz a 

V^y R 


ErAr + ZrE^Arp + Erioz 


R" 


R" 


R^ 


> (E15) 


N 


G 24 - ^ 5 (x - jZ)fErIoj 2 ^ - ^rr?ErIoxz 

j.l V 


R" 


R^ 


N 

G 33 = 2 
j=l 

N 

G34 = ^ 6(x - ji) 

j=l 

N 

G44 = ^ 5(x - ji) 

j=l 


— + 2 Zj»EpAp — - + Ej«Io 2 ;| — j 


®rIoxz ^ " ^r^ErlQx^fi:] 


E„I 


r X 


+ z 


R‘ 


- 2 zrr 7 E],I- + GfJ; 


2 r ' ^ xlR 
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APPENDIX E - Continued 

where 

77 = 1 - ^ ^Note that Ej.Io2; = Eylr + 

In this form, the elements of the ring stiffness matrix G can be compared with results 
obtained in other buckling studies of ring-stiffened shells. (See refs. 28 to 30.) 

Equilibrium Equations 

The condition for vanishing of the first variation of the ring energy (oiIrB = o) yields 
the contribution of the more accurate ring theory to the buckling equilibrium equations. 

If the terms in 6U, 6V, and 6W are collected, the set of buckling equations which are 
the counterpart of equations (C15) with a more accurate ring theory can be written as; 

AllU" + ^Ai 2 - Gii)u + A13V - G12V + (ai 4 - TaAis^W” + (Ai 6 - Gi4 - TAA^jw 

(Ai8 " ^13 " ’^aAi9)w + AiioM' = ^ (E16a) 

A21U' - G12U + A22V" + (a 23 - P2i)v + A24W" + (A25 - G24 - TaA 26)W' 

+ (a 27 - P22 - TaA 28 )w = 0 (E16b) 

(A3I - Gi 4 - TaA 32 )u’ + r3iU + A33V” - G24V’ + (A34 + P32 - TaA3s)v 
+ (^36 - ^33 - TaA 37 )w” + (A38 - T34 - TAA39 - t| A gi^jw 
+ (^311 - G34 - TaA312)w - M" = 0 (E16c) 

A4iU' + A42V + A43W" + (A44 - TAA45jw’ + A4gW + M = 0 (E16d) 

where 

N 2 

^21 ” ^22 " ^ 

j=l ^ 

N 

^22 = ^23 - X 
j=l 


®rAr ^ + ZyErAr^ 
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APPENDIX E - Concluded 


^31 " ^13 " ^14 



Equations (E16) were cast in finite-difference form using the techniques described 
in appendix C and reference 17; the results were programed as part of the determinant 
plotting solution in the computer program entitled BAMSOC I (see appendix D). However, 
numerical results employing the more accurate ring theory are not available. 
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TABLE I.- STRUCTURAL PROPERTIES FOR BUCKLING CALCULATIONS 


(a) SI Units 


Type of 

Esk 

Msk 

^k 

Er 

Ar 

Ir 

Zr 

GfJr 

Es 

^s 

Is 

zs 

GsJs 

cylinder 

GN/m2 

Per K 

GN/m2 

cm2 

cm4 

cm 

N-m2 

GN/m2 

cm2 

cm4 

cm 

N-m2 

Ring- stiffened 

183 

0.3 

1.43 X 10-5 

183 

0.626 

1.228 

-1.364 

0 






Stringer -stiffened 

72.4 

.32 

2.52 






72.4 

2.394 

3.663 

-2.344 

68.44 

Ring- and 

100 

.30 

0.9 

113 

1.587 

5.732 

-5.060 

9.999 

113 

.319 

.1699 

-.864 

0 

stringer -stiffened 















(b) U.S, Customary Units 


Type of 

Esk 

^sk 


Er 

Ar 

Ir 

*_ 

Zr 

GrJr 

Es 


Is 

♦- 

Zs 

GsJs 

cylinder 

psi 

Per °F 

psi 

in2 

in^ 

in. 

lb-in2 

psi 

in2 

in^ 

in. 

lb-in2 

Ring-stiffened 

26.5 X 10® 

0.3 

7.95 X 10-6 

26.5 X 106 

0.097 

0.0295 

-0.537 

0 






Stringer -stiffened 

10.5 

.32 

14 






10.5 X 106 

0.371 

0.088 

-0.923 

23 800 

Ring- and 

14.5 

.30 

5 

16.4 

.246 

.1377 

-1.992 

3477 

16.4 

.0495 

.004082 

-.342 

0 

stringer -stiff ened 















* Based on reference surface at skin middle surface. 
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RING DETAIL 

Figure 2.- Geometry of ring-stiffened cylindrical shells from reference 11. 

Dimensions are in inches (cm). 
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Figure 3.- Thermal buckling predictions for ring-stiffened cylinders. 
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Figure 5.- Geometry of stringer -stiffened cylinder from references 23 and 24. 

Dimensions are in inches (cm). 
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(a) Longitudinal skin stress. 


Figure 6.- Stress distribution in clamped stringer-stiffened cylinder under 
combinations of axial compression and thermal loadings. (Negative sign 
on stresses denotes compression.) 
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RADIAL DISTANCE, z/h 


00 

o 





I — I 


-150 0 


N^, Ib/in. T F 
(MN/m) 

3725 

(0.652) 0 


-40 0 



5275 

(0.924) 



0 150 

a^,/VlN/m 


50 

(28 '^K) 


BO 

(72 °K) 


(c) stringer stresses. 


Figure 6.- Concluded. 





INTERNAL RINGS 
AND STRINGERS- 



Figure 9.- Geometry of ring- and stringer-stiffened cylinder. Dimensions are in inches (cm). 
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(a) Skin stresses. 

Figure 10.- Stress distribution in simply supported ring- and stringer -stiffened 
cylinder under a combination of axial compression and temperature 
^Nx = 1650 Ib/in., 0.289 MN/m; Tgk = 360° F, 200 K^. 
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(b) Stiffener stresses. 


Figure 10.- Concluded. 
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AXIAL LOAD APPLIED AT SKIN MID-PLANE 



Figure 11.- Buckling interaction curve for simply supported ring- and stringer-stiffened cylinder under 
combinations of axial compressive, thermal, and pressure loadings. 
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Figure 12.- Radial buckling deformation in thermally loaded ring- and stringer -stiffened cylinder. 
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